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ABSTRACT 

The   Provincial    Water   Management   Goals  and   Objectives   provide   for  a 
mixing  zone   and  a  zone  of   passage    in   the  vicinity  of  an  outfall    in   a 
receiving   water   body  to   accommodate  the   conflicting   needs  of 
effl  jent  discharge   and   habitat    requirements  of  aquatic   biota.      This 
report  deals  with  descriptions   of  field   survey  and   prediction 
procedures   applicable   to  mixing   zones   below   bank   outfalls    in    shallow 
rivers.      It   is   mainly  an   abridged   version   of  the  MOE  publication, 
Water  Resources   Paper  #14. 

The  mathematical   model    formulations,  based  on  the  stream  tube 
concepc    (which   considers   partial    cumulative  discharge   instead  of  the 
lateral    distance  as  the  transverse  co-ordinate),   are  applicable  to  a 
natural    stream   stretch   in  which   the  discharge   remains   constant.      The 
model    is   developed   from  an   analytical    solution  of  the   steady   state, 
2-D  convection-diffusion  equation,   modified   to   account   for  the 
longitudinal    variabilities   of  decay  and  dispersion   parameters. 

The  concentration  of  a  pollutant   along   a  given   lateral    boundary  of 
mixing   zone    (or   limited   use   zone,   LUZ)    reaches   a  maximum  at   some 
longitudinal    distance   below  the  outfall;   this   point    is   termed 
"critical    point".      Analytical    expressions   are   presented    for  the 
co-ordinates  of  the   critical    point,    as  well    as   the   allowable 
effluent   concentration   to  meet   a   specified   instream  concentration 
criterion.     A   family  of  curves    is   presented   for  a    range  of 
dimensionless   co-ordinates   of  the  critical    point. 

Procedures   for   carrying  out    field    studies,   data   analysis   and 
parameter   estimation   are   presented.      A  step-by-step   design 
procedure,    based   on   the   critical    point  method,       is  outlined   to 
compute   critical    concentration   and   allowable   effluent   concentration 
of  a   pollutant,    as  well    as   the   longitudinal    boundary  of  a   LUZ.      An 
example    illustrates   the   application   of  the  design   procedure.     The 
computations   can   be   performed   through   the   use  of  charts   and  desk-top 
calculators.     The  application   of  model    to   assess   the   impact   of 
typical   management  options   on   allowable   effluent   concentrations   of 
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total  anmonia  and  associated  maximum  longitudinal  spreads  is  also 
illustrated.  The  results  indicate  that,  in  general,  the  allowable 
effluent  concentration  as  well  as  the  maximum  longitudinal  spread  of 
a  LUZ  would  increase  due  to:   (a)  an  increase  in  the  lateral 
boundary  of  the  LUZ;  (b)  an  increase  in  the  streamflow  rate;  and  (c) 
a  decrease  in  temper-ature. 

Data  on  the  distribution  of  residual  chlorine  in  tne  Soyne  River 
below  Alliston  were  used  to  validate  tne  model  satisfactorily  for  a 
nonconservati ve  material.  For  comparison  purposes,  allowable 
effluent  concentrations  of  residual  chlorine  to  meet  a  specified 
instream  criterion  were  predicted  by  the  LUZ  concept  and  on  the 
assumption  of  instantaneous  complete  mixing  (ICM)  of  effluent  with 
the  streamwater.  The  results  indicate  that  in  some  cases,  the  ICM 
assumption  or  the  dilution  ratio  concept  may  lead  to  a  LUZ  with  too 
large  a  lateral  boundary  (and  consequently,  too  small  a  zone  of 
passage)  due  to  an  underestimation  of  treatment  required  in 
comparison  to  the  predictions  based  on  the  LUZ  concept;  whereas,  in 
some  other  situations,  the  predictions  based  on  the  former  method 
may  result  in  treatment  requirements  that  are  too  stringent  in 
comparison  to  those  of  the  latter  method.  Thus,  in  the  case  of  a 
mixing  zone  below  a  bank  outfall,  in  which  lateral  mixing  of 
effluent  is  gradual,  it  is  inappropriate  to  determine  the  allowable 
effluent  concentration  of  a  pollutant  based  on  the  ICM  assumption. 


Ill  - 


ACKNOWLEDGEMENTS 

The  author  is  thankful  to  the  Chairman,  Secretary  and  Members  of 
Working  Group  II,  Water  Management  in  Ontario-Implementation,  and  to 
the  staff  members  of  the  Water  Resources  Branch  and  Regional  Offices 
of  the  Ministry  of  the  Environment,  who  provided  constructive 
comments  and  suggestions  during  various  stages  of  the  preparation  of 
this  report. 


TV  - 


TABLE  OF  CONTENTS 

Page 

ABSTRACT  i  i 

ACKNOWLEDGEMENTS  iv 

TABLE  OF  CONTENTS  v 

LIST  OF  FIGURES  vii 

LIST  OF  TABLES  vn 

LIST  OF  NOTATIONS  viii 

I.  INTRODUCTION 

1.1  General    1 

1.2  Mixing   Zones   Below  Bank   Outfalls-Definitions    3 

1.3  Lateral    Boundaries  of  LUZ  and   ZOP   5 

1.4  Objectives    5 

II.  MODELLING  CONCEPTS 

2.1  Description   of  Model    8 

2.1.1  Assumptions    8 

2.1.2  Mathematical    Formulations    9 

2.1.3  Graphical    Solution    14 

2.1.4  Computer   Program  MIXCALBN 16 

2.2  Limited   Use   Zone   Boundary   and   Allowable  Effluent 
Concentration    17 

2.3  Critical    Point  Method    19 

2.3.1  Critical    Co-ordinates    19 

2.3.2  Graphical    Method    20 

2.3.3  Effects   of  Neglecting   Side-wall 

Reflection   of  Material    20 

2.3.4  Maximum  Longitudinal    Spread    22 

2.4  Effects  of   Background  Concentrations    25 

2.5  Scale-up  of  Model    Parameters    26 

2.6  Computer   Program  MI XAPPLN 27 

III.  DATA  COLLECTION  AND   ANALYSIS   PROCEDURES 

3.1     Field   Survey   Procedure    29 


-  V   - 


TABLE  OF  CONTENTS  (continued) 

Page 

3.2  Data   Analysis   Procedure    30 

3.2.1  Mean   Depth   and   Velocity    31 

3.2.2  Transverse  Distribution  of   Flow 32 

3.2.3  Mass   Flux   and   Relative  Concentrations    . 34 

3.3  Simplified  Method   for  Estimation  of    i  35 

3.4  Computer   Program  for   Data   Analysis    36 

IV.  APPLICATION  OF   PREDICTION   TECHNIQUES 

4.1  General    Procedure    37 

4.2  Design   Procedure   for   the  Critical    Point  Method    38 

4.3  Total    Aimonia   in   the  Mixing   Zone  of  Grand   River 

below  Waterloo    39 

4.3.1  Illustration   of   Step-by-Step   Design 

Procedure    40 

4.3.2  Management   Options    43 

4.3.3  Modell ing  Results    44 

4.4  Residual    Chlorine   in   the   Boyne   River  below 

All  1st  on 47 

4.4.1  General    Description   of   Studies    47 

4.4.2  Data  Analysis   and   Parameter  Estimation    50 

4.4.3  Model    Validation    Study    52 

4.4.4  Critical    Point  Method   Predictions    53 

V.  SUNWARY  AND  CONCLUSIONS    56 

REFERENCES    58 


-    VI     - 


LIST  OF   FIGURES 


Page 


Fig.      1     Configurations  of  Mixing   Zone   Boundaries    in   Rivers    ..  2 

Fig.      2     Schematic   View  of  Mixing   Zone    4 

Fig.      3     Effect  of   Cross-Sectional    Flow  Distributions   on 

Lateral    Boundary  of   LUZ  and   ZOP    6 

Fig.     4     Dimensionless   Concentration   Distribution 

in   Mixing   Zone    15 

Fig.      5     Critical    Point  of   LUZ    18 

Fig.      6     Non -Dimensional    Co-ordinates   of  Critical    Point 21 

Fig.      7     Side-wall    Reflection   Effects    23 

Fig.     8     Graphical    Method    for   Longitudinal    Spread   of   LUZ   24 

Fig.     9     Determination   of  Longitudinal    Spread-Design   Problem   .  42 

Fig.    10     Schematic    Layout   of   Study  Area   -   Boyne   River    48 

Fig.    11     Plot   of    (c,/c       )      versus    (x/b)    -  Boyne   River    .  51 

Fig.    12     Attenuation   and   Decay  of   TRC    in   the   Boyne   River    51 

Fig.    13     Model    Validation   Results   -  Boyne   River    54 

Fig.    14.    Distribution  of   TRC  along   the   Lateral 

Boundary  of   LUZ   54 


LIST  OF   TABLES 


Page 


Table  1     Surmiary  of  Nondimensional    Coefficients,    b  ,    in 

Natural    Streams   U 

Table  2     Shoreline  Concentration   Distribution   - 

Design    Problem   41 

Table   3     Arnnoma    in   the  Grand   River   -  Sumnary  of  Predictions    .  45 

Table   4     Surrmary  of  Boyne   River   Field   Survey  Data   - 

August    11-12,    1976    49 

Table   5     Results   of   Critical    Point   Method    53 


-    VI  1 


LIST  OF   NOTATIONS 
The   following   notations   are   used    in   this   report: 

A:  Cross-sectional    area   of  stream  channel  . 

B:  Weighted  mean  value  of  width  of   stream  channel. 

b:  Top  width   of  water   surface   at   a  given   transect. 

bex:  Width   exponent    in   Leopold-Maddock   Equation. 

c:  Concentration   of  a  material    at   a   given   point   in  mixing   zone, 

c    :  Completely  mixed   instream  concentration   of  a  material. 

a 

c.  :  Background  concentration  in  river  water. 

c  :  Concentration  of  a  conservative  material  at  a  given  point. 

c  :  Effluent  concentration  of  a  material. 
e 

c.  :  Depth-averaged  concentration  at  a  lateral  distance,  y,  , 

from  outfall  shoreline. 

c.  :  Critical  concentration. 

C|^:  Ratio  (c^/c^). 

c„:  Ratio  (c,  /c^  ). 

r  La 

c    :  Specified    instream  concentration   criterion. 

c.  :  Concentration  of  total    aimonia. 

c   :  Concentration   of  un-ionized   amnonia 

c    .:  Allowable   effluent   concentration 
eA 

c^^:  Ratio   (c/Catr) 

c  ,:  Instream  concentration  at  a  downstream  boundary  due  to 
wa  •' 

effl uent  discharge, 

c^  :  Average  concentration  of  material  at  a  given  transect 

obtained  from  (mass  flux/discharge). 

c   :  Maximum  concentration  of  material  at  a  given  transect, 
c   .  :   Critical  concentration  estimated  during  iteration  procedure. 

D  :  Diffusion  factor. 

E. :  Nondimensional  term  defined  in  Eq.  11. 

F:  Net  mass  flux  at  a  given  transect. 

fj^:  Mass  flux  of  material  in  a  vertical  strip,  k. 

f,:  Constant  in  the  resistance  equation  (Eq.  25). 

^2'-  Exponent  in  Eq.  25. 

G,  :  Nondimensional  term  defined  in  Eq.  11. 

h:  Local  depth  of  water  at  a  lateral  distance,  y. 

H:  Mean  depth  in  a  given  stretch  of  river. 
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LIST  OF   NOTATIONS    (continued) 


h:  Mean   depth   at   a  given   transect. 

h,  :  Local    depth   at   a  distance,  y,  ,    from  outfall    shoreline. 

hex:  Depth   exponent    in   Leopold-Maddock   Equation. 

ka:  Dissociation   constant   for  armonia. 

k,  :  Decay   rate  associated   with  the   background   concentration   of 

nonconservati ve  material. 

k  ,  :  Decay   rate  of  nonconservati ve  material    in    river. 

k    :  Attenuation    rate  of  a  material    due  to   transverse   dispersion. 
P 

k    :  Gross  decay   rate    (due  to   dispersion   and   chemical    reaction). 

k      :  Weighted  mean   value  of  the  decay   rate,   k,. 
av  ^  d 

m^ :  Average  value  of  metric   coefficient    (or   scaling   factor) 
associated  with  x-axis. 

n:  Number   of   images. 

n':  Manning's    roughness   coefficient. 

p:  Ratio    (q/Q). 

p^_:  Ratio    (q^/Q). 

p    :  Maximum  value  of  p  at  which   n   is   zero. 
m 

p^:  Ratio    (q3/Q). 

pH:  pH  val ue  of  water. 

pka:  log    (1/ka). 

Q:  Discharge   in   river   below  outfall. 

Q   :  Effluent   flow  rate. 

q:  Cumulative  partial    discharge  measured   from  a   reference   bank 

of   stream  channel . 

q. :  Lateral    boundary  of  limited   use   zone. 

q    :  Cumulative  partial    discharge  between   reference   bank   and 

effluent  discharge  point    in    river. 

R    :  Factor  derived    from  a  chemical    equilibrium  relationship. 

R.:  Decay   factor    in    i^h   ^^^^^  Qf  stream. 

s    :  Energy  gradient. 

T°:  Temperature   of  stream  water. 

t:  Time  of   travel  . 

U:  Mean   velocity   of   flow   in   a   stretch  of  river. 
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LIST  OF   NOTATIONS    (continued; 


u:  Depth-averaged    local    velocity   at   a   lateral    distance,  y. 

u.  :  Depth-averaged   local    velocity   at   a   distance,  y,  ,    from 

outfall    shorel ine. 
a:  Mean   velocity   at   a   given   transect. 

uex:  Velocity  exponent    in   the  Leopold-Maddock   equation. 

x:  Longitudinal    distance  below  outfall. 

X.  :  Critical    distance. 

X   :  Longitudinal    distance   from  outfall    to   a  downstream  boundary. 

X        :         Longitudinal    distance  along  the  discharge  shoreline  where 
see 

the  concentration    is   c    .    as   a   result  of  a  given   effluent 

concentration,   c   . 
e 

X      ■ ^  :       Critical    distance  estimated  during   iteration   procedure, 
cnt  ^ 

X        .:        Value  of  x         when   effluent   concentration   is   c    ,. 
sceA  see  eA 

y:  Lateral    distance  measured    from  outfall    shoreline. 

S  :  Dimensionless  coefficient. 

(p  :  Dimensionless   term   (B   x/B). 

Ç,:  Ratio   (G|_/E|_). 

\()  :  Shape-velocity    factor. 

0:  Temperature  correction   factor. 


I.      INTRODUCTION 

1.1      General 

Wastewater  effluents  discharged  into  streams  and  rivers  contain 
pollutants  such  as  residual  chlorine,  amnonia  and  phenol  which  have 
the  potential  to  be  toxic  to  aquatic  biota.  The  toxic  effects  of 
such  pollutants  are  minirmzed  by  suitable  siting  of  the  outfalls  in 
receiving  water  bodies.  The  types  of  outfalls  that  may  be  employed 
in  shallow  streams  and  rivers  are  described  briefly  in  the  following 
paragraphs.   (Note:  In  this  report,  streams  and  rivers  with  width  to 
depth  ratios  of  10  or  more,  are  considered  to  be  shallow.) 

(a)  A  diffuser  outfall  extending  across  the  entire  width  of 
river  channel  so  as  to  achieve  instantaneous  complete 
mixing  of  effluent  with  river  water  (see  Fig.  la).   In  this 
case,  the  allowable  effluent  concentration  of  a  pollutant 
Is  determined  from  simple  dilution  computation. 

(b)  Pipe  outfall  located  at  a  bank  as  in  Fig.  lb,  so  as  to 
create  an  effluent  plume  that  hugs  the  shoreline  for  some 
distance  downstream  of  discharge.  Thus,  a  portion  of  the 
cross-section  of  the  stream  is  maintained  to  serve  as 
habitat  for  fish  and  other  biota.  This  type  of  outfall  is 
usually  known  as  a  bank  outfall,  and  is  commonly  employed 
in  shallow  streams.   In  this  case,  the  allowable  effluent 
concentration  is  dependent  on  the  transverse  dispersion 
characteristics  of  the  river  channel  and  decay  of  the 
pollutant  in  the  nver.  Assessment  technique  for  this  case 
is  somewhat  more  complex  than  the  previous  one;  and  it  is 
the  subject  matter  of  this  report. 

(c)  Pipe  outfalls  located  in  a  cross-section  of  "iver  at  some 
point  other  than  a  bank  (see  Fig.  Ic).   In  this  case,  the 
effluent  plume  would  be  somewhat  wider  and  the  longitudinal 
extent  would  be  smaller  than  that  of  case  (b)  due  to  the 
spreading  of  effluent  on  either  side  of  the  centre  of  plume. 
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FIGURE  1 ,   CONFIGURATIONS  OF  MIXING  ZONE  BOUNDARIES 
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(d)  Diffuser  outfalls  (see  Figs.  Id  and  le):  The  initial 

dilution  is  higher  in  these  cases  in  comparison  to  the  pipe 
outfall  cases  described  above.  The  plume  would  be  wider, 
but  of  shorter  length.  Assessment  techniques  for  these 
cases  are  also  complex,  and  are  somewhat  similar  to  those 
for  the  pipe  outfall  cases.  However,  when  the  diffuser 
outfall  extends  right  across  the  river  width,  instantaneous 
complete  mixing  is  achieved  (see  case  (a),  described  above). 

1.2     Mixing  Zones  Below  Bank  Outfalls  -  Definitions 

When  effluents  are  discharged  through  bank  outfalls,  the  plume  tends 
to  hug  the  shoreline  for  some  distance.  Vertical  mixing  is  usually 
complete  within  a  short  distance  below  the  outfall  (about  50-100 
times  the  channel  depth).  The  effluent  spreading  is  gradual, 
resulting  in  concentration  gradients  being  set  up  in  the  lateral 
direction.  Ultimately,  the  cross-sectional  concentration 
distribution  attains  uniformity  at  some  distance  below  the  outfall. 
The  zone  of  the  receiving  stream  between  the  outfall  and  the  nearest 
cross-section  at  which  the  concentration  distribution  becomes 
uniform  is  generally  known  as  "mixing  zone",  and  the  corresponding 
distance  is  termed  "mixing  zone  length".  The  distance  below  the 
outfall  at  which  the  effluent  spreads  across  the  entire  width  of  the 
river  is  termed  "crossing  distance".  Within  the  mixing  zone,  a 
portion  of  the  cross-section  of  the  river  in  which  the  concentration 
of  a  pollutant  of  concern  may  not  comply  with  a  specified  instream 
criterion  is  termed  the  "limited  use  zone"  (LUZ);  the  remaining 
portion,  which  is  designated  to  serve  as  suitable  habitat  for  fish 
and  other  desirable  aquatic  life  is  known  as  the  "zone  of  passage" 
(ZOP)  (9,  10)^ 

The  salient  features  of  a  mixing  zone  are   depicted  in  Fig.  2. 
Usually,  the  location  of  the  lateral  boundary  is  designated  as  a 
fraction  of  channel  width  or  river  discharge  for  simplicity,  which 
indicates  that  the  shape  of  LUZ  would  be  a  rectangle.   However,  the 
concentration  decreases  in  the  downstream  direction  due  to 
dispersion  and  decay,  and  hence,  the  lateral  boundary  of  LUZ  will  he 
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FIGURE   2.      SCHEMATIC  VIEW  OF  MIXING  ZONE 
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located  along  a  specified  criterion  concentration  contour  line  as 
indicated  by  the  broken  line  in  Fig.  2. 

A  limited  use  zone,  as  defined  above,  denotes  a  zone  of  a  receiving 

water  body  that  has  traditionally  been  known  as  a  mixing  zone  in 

water  quality  practice.  The  broader  definition  of  a  mixing  zona, 

given  in  the  foregoing  paragraph  (incorporating  LUZ  and  ZOP),  is 
adopted  in  this  report. 
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Lateral    Boundaries  of  LUZ   and   ZOP 


In  the  case   of  a   rectangular  channel  ,   the   lateral    boundary  of   LUZ 
can   be  designated   to   be  a   fraction   of  total    channel    width   since   tne 
discharge  per   unit  width   of  channel    is   the   same.      However,    in   a 
natural    river,    the  channel    widths,   depth   profiles   and   velocity 
distributions   at   successive  cross-sections   are   generally  nonuniform 
as   shown    in   Fig.    3;    in   some   locations   stagnant   areas  may   be 
present.     Consequently,    the  discharge  per   unit   width  of  channel    at 
various   cross-sections  will    vary,    as   illustrated    in   Fig.    3  by   the 
cross-sectional    flow  distribution   curves.     For  example,   40%  of 
discharge   occurs    in  widths  of   55%,    35%  and   30%  at   the   transects   A,    B 
and   C,    respectively;   on   the  other   hand,  within   a   40%  width,   the   flow 
rates   are   16%,    45%  and   65%  at   the   three   successive   transects.      This 
example    indicates   that    it    is  more   appropriate   to  designate  the 
lateral    boundary  by  a   fraction   of   total    river   flow   (rather   than   by   a 
fraction   of  total    width).      Such   a  designation   of  the   lateral 
boundary  will    account   for   the   effect   of  variations    in   channel 
geometrical    and   hydraulic   parameters   at   different   cross-sections   in 
a  mixing   zone. 


1.4 


Objectives 


The  mam  goal  of  this  report  is  to  present  field  procedures  and 
prediction  techniques  applicable  to  assess  water  quality  in  mixing 
zones  of  shallow  rivers  and  streams  below  bank  outfalls.  The 
specific  objectives  of  the  study  are  as  follows: 
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(i)      To   present   a   summary  of  the  modelling   concepts; 
(ii)      To  describe    field   data   collection   procedures; 
(iii)      To  outline  data   analysis   and   parameter   estimation 

procedures  ; 
(iv)     To  present   a   step-by-step   design   procedure   including   an 
i 1 1 ustrati  ve   example; 
(v)     To   illustrate   the   impact   of  some  management   options   on   the 
instream  concentrations,    allowable   effluent   concentrations, 
and    longitudinal    and   lateral    boundaries   of  mixing   zone 
(with   respect   to  ammonia   and   total    residual    chlorine). 
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II.  MODELLING  CONCEPTS 

2.1  Description  of  Model 

The  mathematical    model    is   based   on   the   "stream  tube"   concept.     This 
model,   developed  by   Yotsukura   and   Cobb    (16),    considers   the 
cumulative  partial    discharge,   q,    at   a   given   cross-section    instead  of 
the   lateral    distance,   y,    as    independent   variable.      In   this   approach, 
the   cross   section    is   divided   into   a  number   of  vertical    strips   termed 
"stream   tubes",    such   that   the  discharge  within   each   stream   tube    is 
the   same;  whereas,    in   the  traditional    approach,   the   strips   are   of 
equal   width.      Thus,    the   cross    sectional    concentration  distributions, 
c    (x,   q),   predicted   by  the   stream  tube  model    will    be   functions   of 
q.      These  distributions   can   be   transformed    into   the   traditional 
c(x,  y)   versus  y  plots   by  knowing   the   relation   between   q  and  y  at 
each   transect.      The  mathematical    formulations  of  the   stream  tube 
model    have  been   presented   by   Yotsukura.  and   his   co-workers    (16,    17), 
and   the   derivation   of   critical    point  method   from   those   formulations 
has   been   presented    in   the  Water   Resources    Paper   #14   (Ref.   8).     A 
sunmary  of   those  developments    is   presented    in   this   chapter. 

2.1.1  Assumptions 

The  major   assumptions   of  the  model    are   as   follows: 

1.  The  density  of  effluent    (or   solute)    is   the   same   as   that  of 
the   receiving  water.      This   assumption    is   reasonably 
satisfactory   for  most  of  the  municipal    effluent  discharges 
to   surface  water   bodies. 

2.  The   concentration   distributions    in   the    far-field   are   not 
affected  by    the   near-field  mixing   processes    (eg.    dilution 
due   to    initial    momentum  of  jet).      Usually,    the   jet-induced 
diffusion   approaches   the  ambient   diffusion    in   a    short 
distance   below  a   source    in   a   shallow  river. 


3.  The   vertical    distribution   of  effluent    (or   solute)    in   the 
river   channel    is   uniform.      Generally,    the    longitudinal 
distance    required   to   attain   vertical    uniformity    is   fairly 
short    in    shallow    rivers,    being   on   the   order   of   50   to    100 
times   the   channel    depth. 

4.  The   transport  due   to    longitudinal    dispersion    is 
negligible.      In   the   case   of  continuous  effluent   discharge, 
this   transport    is   very   small    in    relation   to   that  due   to 
convection   and    lateral    dispersion,   thus   justifying   the 
assumption. 


2.1.2 


Mathematical    Formulations: 


The    simplified  two-dimensional    convection-diffusion  equation, 
describing   the   distribution   of  a   nonconservati  ve  material    in   the 
far-field    region   of   the  mixing  zone    in   a    shallow   river   under   steady 
state   conditions,   can   be   written    in   the   following    form    (assuming 
that   the   decay  of  the  material    follows      first-order   kinetics): 


3x 


D 


;52 
o  c 


k  ,c 
d 


(1) 


in  which   x   =  longitudinal    distance   below  the   source;   q   =  cumulative 
partial    discharge  measured   from  a   reference    bank;   c   =  concentration 
of  pollutant   at   a   point    (x,   q);   U   =  average   velocity   of   flow   in   the 
x-direction;   k.    =  first-order  decay   rate   coefficient    (assumed   to 
be    independent   of  x  and  y);   and   D     =  diffusion   factor. 


The   diffusion    factor    is   assumed   to   be    constant    at   a   cross-section; 

but,    it    is   dependent   on   the    variations    in   the   width   and   alignment 

of    river  channel,    transverse   dispersion   coefficient   and 

cross-sectional    distributions   of  depths   and   velocities.      In   a 

2 
previous    publication    (8),    it    is   shown  that    D     =    g  Q  /B   and    B  = 

m     ,|;    g        ,   where   Q   =   river   discharge   just   below  the   outfall;      B   = 
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average  channel    width;     6=  nondimensional    coefficient    (or 
nondimensional    diffusion   factor);   m     -  average  value  of  metric 
coefficient    (or  scaling    factor),   which   accounts   for  change   in  width 
and   alignment   of  channel    in   the   x-direction;     4)=  shape-velocity 
factor,   dependent   on   the  deviation  of  cross-sectional    mean   depth   and 
velocity   values    from  their   local    values;   and   B      =  nondimensional 
transverse  dispersion   coefficient. 


The   parameter  s  (or   D   )   may  be   estimated   from  data   on   tracer 
concentration  distributions   at    successive  transects,    and   channel 
geometric   and   hydraulic   data    (1,    8,    16);   a   simplified  method   of 
estimating   8    is  given    in   the   next   chapter.      In   the  absence  of   such 
data,   an   approximate   value   for   a  given   river   stretch  may  be   selected 
from   results  of   studies   carried  out    in    rivers   and   streams 
elsewhere.      Table    1  shows   a   sunmary  of  the   parameter   values, 
obtained   from  data   published    in   the   literature    (1,   8,    16).      From 
this   table,   the   values   of     S   ^.r^   seen   to   be    in   the   range   10' 
10     .      Generally,    lower   values   are  observed    in   slow-moving 
straight  channels  with   low   turbulence   levels,   whereas   higher  values 
^r^  measured    in   highly  turbulent,  meandering   streams. 

The  concentration  of  a   pollutant   at  a  given   point    in   the  mixing   zone 
may  be  obtained   from  a  numerical    solution    (eg.    finite-difference 
method)    or   from   hn   analytical    solution   of  Eq.    1,      The  numerical 
solution   techniques   require  detailed  data   on   channel    morphology  and 
velocity  distributions,    and   the   calculations   àr^   performed  with   the 
aid  of  a  digital    high  speed   computer;   therefore,   numerical    solutions 
àVQ  employed  mostly    in   complex    situations.      In   contrast,    analytical 
solutions   require   lesser  amount   of  data,   and   computations    involved 
à?Q  much    simpler.      When    the   longitudinal    variabilities   of   geometric 
and   hydraulic   parameters   of  channel    and   pollutant  decay   rates    in   a 
given    river    stretch   are   significant,    the   predictive   capabilities  of 
the   analytical    solutions   can   be    improved   through  the   use   of  moving 
average   values   of   the   pertinent    parameters   as    shown    in   the  Water 
Resources    Paper   #14   (Ref.    8).      The   computer   programs   presented    in 
the  Water   Resources   Paper   #14   àVQ   based   on   the   use  of  moving   average 
values   of   various   parameters.      Since  manual    computation   of  moving 
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averages   are   too   laborious,   the   average  parameters   of  the   river 
stretch   under  consideration  may   be   used   as   an   approximation. 
Therefore,   the   analytical    solutions   presented   in   this   report   are 
written    in   terms  of  average  values  of  the   appropriate   parameters. 

Analytical    solutions   of   Eq.l   applicabU:   to   continuous  discharge  of 
effluents   from  pipe   and  diffuser  outfa   Is,    are  given    in   the 
literature    (16).      For   the  case  of  mixing   zones   below  bank   outfalls 
in   shallow   streams   and    rivers,   a   particular   solution  of   Eq.l   can   be 
written    in   the   following   dimensionl  ess   form   (8,    13): 

+  00 

Cr   (  *.,P)    =  R^    (   ^^  S^''^     Y,   exp|-(p-2n)^/4   ^\,  (2) 

n=-<»        ^  ^ 

in  which  the  dimensionless  terms  are  defined  by 

cp.  =  S.x./B^   ;   p  =  q/Q; 

Cp  (  ^  ^,  p)  =  c(x. ,  q)/c^;    c^  =  c^Q^/Q; 

R.    =  exp   (-k.    x./U.)  (3) 

The  variables   appearing    in   Eqs.2  and   3  are   as   follows: 

X.        =     Longitudinal    distance   from  outfall    to   a   transect,    i. 

D 

i   =  Average  channel  width  of  the  river  stretch  of  length, 

X.  , 

1 
U.   =  Average  velocity  of  flow  in  the  river  stretch. 

q   =  Partial  cumulative  discharge  at  a  cross-section, 

measured  from  the  outfall  shoreline. 

Q   =  Discharge  in  the  river  just  below  the  outfall. 

Q    =  Effluent  discharge. 
e  ^ 

c    =  Concentration  of  pollutant  in  the  effluent. 

g.   =  Nondi mensi onal  coefficient  applicable  to  the  river 

stretch  of  length  ,  x . . 
k.   =  Average  decay  rate  coefficient  in  the  river  stretch  of 

length  ,  x- . 
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c(x. ,q; 


R. 


=  Concentration  of  pollutant  at  a  point  (x- ,q)  with 
respect  to  the  point  of  effluent  discharge. 

=  Average  concentration  in  the  cross-section  of  river 
just  below  the  outfall. 

=  Fraction  of  total  river  discharge. 

=  Di  mensi  on!  ess  concentration,  c(x.  ,q)/c;  . 

1       d 

=  Decay  factor  applicable  to  the  river  stretch  of 

length  ,  x. . 
=  Number  of  images  required  to  account  for  the 

reflection  of  material  from  the  channel  banks. 


By  following  an  empirical  procedure  described  in  previous 
publications  (6,  7,  8),  the  following  expression  for  n   is  obtained: 


n  =  0.5p  +  a.OSVï'. 


(4) 


In   computing   the   concentration    at   a   given   point    (x,q)    from   Eq.2,    the 
sunmation  of  the  exponential    series   terms   is   carried   out    from     n   =-n, 
to   n   =  +n„,   where   n,    and   n„  are   obtained   from  Eq.4  with   the 
negative   and   positive    signs,    respectively. 

Equation   2   is   applicable   to   conservative  materials   (when   k,    is   set 

equal    to  zero),    and   to   nonconservati ve  pollutants  whose  decay 

follows   the   first-order   kinetics.      However,   the   concentrations   of 

some  pollutants   (e.g.   un-ionized  ammonia)    are   also  dependent   on   the 

chemical    equilibrium  of   solutions.      In   such   cases,   the   pollutant 

concentrations    are   obtained   by  multiplying   the   predictions  of   Eq.    2 

by  yet   another   factor,   R    ,    derived    from  chemical    equilibrium 

considerations.      For   example,    in    the   case   of  armionia,    the 

predictions   of   Eq.    2  vvould    represent   the   concentrations   of  total 

aiimonia    (i.e.    the   sum  of   ionized   and   un-iomzed   forms   of  ammonia). 

The   concentration   of  un-ionized    armionia,   c^j,    is   a    function   of   pH, 

temperature,   T°C,    and   concentration   of  total    armionia,   c.  .      The 

relation   between   c     and   c.    can   be   written    in   the    form  c     = 

u  t  u 

R  c    ,   where   the    factor   R      is   given   by   (6,    7,    14): 
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R     =   1.0/(1   +   IQP^^   "  P"),  (5) 


where  pka   =  0.09018  +   (2729.92/(T°  +  273.2)), 


in   which  ka   =  dissociation   constant   for  ammonia,    and   pka   =  -log   kd. 
This   development    is   based   on   the   assumption   that   the   instream 
temperature   and   pH  distributions   at   various   cross-sections   are 
uniform.      This   assumption    is    reasonably   satisfactory  for  most 
practical    situations. 

2.1.3  Graphical    Solution: 

For   a  conservative  material,   the   decay  function   becomes   unity   in 
Eq.    2   (since   k  ,  =0   in   each   reach);   then    it    is   possible   to   calculate 
the   values   of  Cp   for   given   values   of     ^  and   p.     A  graphical 
solution   of   Eq.    2,    obtained    for   a   range  of  values   of     cp.   normally 
encountered    in   practice  with   p=0  to   1.0,    is   presented    in   Fig.    4. 
This   graph   can   be   used   to   obtain   the   concentration   of  a   conservative 
material    at   any   desired   point    in   the  mixing   zone.      The   corresponding 
concentration   of  a   nonconservati ve  material    subjected   to   first-order 
decay   can   then   be   obtained  by  multiplying   each   concentration   value 
by  the   pertinent   decay   factor.     When   the  widths,   velocities   and 
decay    rates   are   reach-dependent,    the   computations   a^e   laborious   and 
time-consuming.      However,   through  the   use  of  weighted  mean   values 
k      ,   B   and    U   for   the   study   stretch   as   an   approximation,   and   by 
estimating    3   as  outlined   elsewhere    in   this    report,    the   following 
procedure   can   be   used   to   obtain   concentration   distributions    from 
Fig.    4: 


1.        Select   the  design   case   parameters,   viz.,    Q,    Q   ,    c     and 
temperatuf'e. 


Estimate    s   ,   k      ,   B   and   U. 


3.        Compute    f     for   known   x,    and   determine   c,    and   p  using 

a 

appropriate   relationships   given    in   Eq.    3. 
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FIGURE  4.   DIMENSIONLESS  CONCENTRATION  DISTRIBUTION  IN  MIXING  ZONE 
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4.  From  Fig.  4,  read  val  ues,  of  Cp  for  various  values  of  p 
corresponding  to  tp    computed  in  step  3. 

5.  Calculate  the  concentration  of  conservative  material  at 

each  point  (x,q)  from:  c  (x,q)  =  c^  c  . 

c         K  a 

6.  To  determine  the  concentration  of  nonconservati ve  material, 

use  the  equation   c(x,q)  =  c  (x,q)  exp(-k  x/U). 

c  a  V 

The  foregoing  graphical  procedure  should  be  used  in  cases  where 
access  to  more  efficient  computational  procedures  is  limited  and 
when  approximate  estimations  are   acceptable. 

The  derivation  of  Eq.  2  is  subject  to  the  conditions  c  — ►  ^  and 

Q  ►"  0;  thus,  Eq.  2  is  generally  valid  when  Q  is  much 

smaller  than  Q.  When  this  condition  is  not  satisfied,  the  values  of 
c(x,q)  predicec  by  Eq.  2  may  be  higher  than  c  ,  although  such 
situations  do  not  occur  in  real  world.  The  conditions  under  which 
the  predicted  concentrations  may  be  higher  than  the  effluent 
concentration,  c  ,  can  be  examined  with  the  help  of  Fig.  4,  as 
outlined  in  Water  Resources  Paper  #14  (Ref.  8).  Such  anomalous 
predictions  should  be  ignored  in  the  analysis. 

2.1.4    Computer  Program  MIXCALBN: 

A  computer  program  MIXCALBN  (Mixing  zone  -  CALi  BratioN  model ) . 
written  in  FORTRAN  language,  has  been  developed  to  predict 
concentration  distributions  in  mixing  zones  below  bank  outfalls. 
The  program  is  based  on  an  analytical  solution  for  a  point  source 
which  is  given  by  Yotsukura  and  Cobb  (16);  the  predictions  of  the 
Yotsukura-Cobb  equation  are   identical  to  those  of  Eq.  2  for  the 
particular  case  of  bank  outfalls.  The  program  MIXCALBN  is 
documented  in  Water  Resources  Paper  #14  (Ref.  8),  interested  readers 
should  refer  to  that  publication  for  details  of  its  usage.  The 
program  is  particularly  suitable  for  model  validation  studies  (i.e., 
calibration  and  verification  of  model).  A  brief  description  of  the 
program  is  provided  below. 
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The  input  parameters  include:  (i)  longitudinal  distances,  x., below 
the  outfall  to  various  transects  at  which  concentration  predictions 
are  required;  (ii)  top  width,  mean  depth  and  mean  velocity  at  each 
transect,  i;  (iii)  decay  rates  in  successive  reaches  (i.e.,  between 
any  two  consecutive  transects)  in  the  river  stretch;  (iv) 
nondimensional  coefficients,  3.:  (v)  river  discharge  and 
background  concentration  of  pollutant;  and  (vi)  effluent  discharge 
and  associated  pollutant  concentration.  The  output  of  the  model 
includes  lateral  concentration  values  from  q  =  0  to  q  =  Q  at  each 
transect.  The  computations  at  successive  transects  make  use  of 
moving  average  values  of  width,  velocity  and  decay  rate. 

The  parameters  input  to  the  MIXCALBN  program  are  obtained  from  field 
study  and  parameter  estimation  procedures;  these  are  outlined  in  the 
next  chapter. 

2.2  Limited  Use  Zone  Boundary  and  Allowable  Effluent  Concentration 

The  boundary  of  a  LUZ  in  a  shallow  river  is  generally  identified  by 
lateral  and  longitudinal  co-ordinates  with  respect  to  the  outfall. 
Usually,  the  lateral  boundary  of  a  LUZ  is  limited  to  the  range  0.2 
to  0.4  times  the  river  discharge,  Q.  The  cumulative  partial 
discharge  between  the  outfall  bank  and  the  accepted  lateral  boundary 
will  be  denoted  by  the  nondimensional  parameter,  p.=q./Q.  A 
specified  pollutant  concentration  criterion,  c  ,  must  be  met 
within  q. .   For  a  given  set  of  values  of  c  ,  Q  ,  q.  and  Q, 
the  longitudinal  distribution,  c(x. ,  q.),  predicted  by  Eq.  2 
attains  a  maximum  value  at  some  x.  and  then  follows  a  decreasing 
trend  as  shown  in  Fig.  5.   The  point  at  which  the  concentration 
attains  the  maximum  value  is  termed  the  "critical  point".   The 
longitudinal  distance  between  the  outfall  and  the  critical  point  is 
termed  the  "critical  distance",  x.  ,  and  the  maximum  concentration 
is,  termed  the  "critical  concentration",  c,  .   Methods  for  the 
computation  of  x,  and  c,  are   presented  in  Section  2.3. 

Knowing  c.  and  c   the  allowable  effluent  concentration,  c  . 
can  be  calculated  from  the  following  expression: 

c^A  -  ^e   S/^L  (^) 
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Eq.   6   follows   from  the   fact   that   a  change   in   c     in   Eq.    2  would 
result    in    a   proportional    change    in   c,  ,    provided   all    other 
parameters   are   held   constant;    the  derivation  of  Eq.    20  has   been 
presented    in   a   previous   publication    (7).  I 

The  maximum  longitudinal    boundary  of   LUZ,    x       occurs   along   the 

discharge   shoreline    in   the   case  of  bank   outfalls   as   shown    in 

Fig.    2.      Methods  of  determining   the  maximum  spreads   resulting   from 

c     and   c    .   are  outlined    in   the   next   section, 
e  eA 

2.3  Critical    Point  Method 


2.3.1         Critical    Co-Ordinates  : 

The  derivation   of  expressions    for   the   co-ordinates   of  the   critical 
point    is   based   on   the   following   assumptions: 

1.  The  channel    hydraulic   characteristics    (viz.,   width,   depth, 
velocity   and  dispersion   coefficient),    and   the  pollutant   decay 
rate   remain   constant   in   the   study   stretch.     When   the 
longitudinal    variability   of  these   parameters    in   the   study 
stretch  are   significant,   average  values   of  the   parameters   can   be 
utilized    (as    stated   earlier). 

2.  The  effect   of  reflection   of  material    from  the  channel    banks   on 
the   concentration   distribution    is   negligible.      The   limitation 
imposed   by  this   assumption    is   outlined   later   in   this   section. 

The  details  of   derivation   of   the   expression    for   the   critical 
co-ordinates   have  been   described    in   previous   publications    (7,    8). 
For  a    nonconservati ve  pollutant,    the   expressions   for   the   critical 
distance,   x.  ,   and   the   critical    concentration,   c,  ,   are  as    follows: 


U       I  ^  5 '  \ 

-1      +    V   1   +  4    (p-   k.B/U^  )     ] 


X      =  ^^       -1      +    V   1   +  4    (pr   k   B/U^)     ]  (7) 

^        4k       '  ^     ^ 


d 


r 


k.x,  pi^  B 


Cl   -   c^      '    — exp  r-  :i^     -  li 1  (8) 
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For  conservative  materials,  expressions  for  x.  and  c.  can  be 

derived  from  Eq.  2  by  setting  k  ,  =  0  and  neglecting  the  effect  of 

2 
images.  The  resulting  expressions  are  x,  =  B  p./2  s  ,  and 

Cl  =  C3/2.07p^. 

2.3.2         Graphical    Method: 


The  expressions    for   the   co-ordinates   of  the   critical    point   given   by 
Eqs.    7   and   8  can   be  written    in  dimensionl ess   terms   as   follows: 


\ 

0. 

J 

,25  <  -1 

'^r 

TT    G. 

|-1   +V  1   +  4  p^^   E^.    )       , 


(9) 


(10) 


where   the   dimensionless   terms  are  defined   by: 


E^   =  By  SU;         9|_   =  G^/E^  (li; 

Figure   6   shows   plots   of   G,    vs.   E,    and   c     vs.   E,    for   various 
values   of  p,    ranging   from  0.1  to   0.6.      The   family  of  curves   shown 
on   this   graph   can   be   used   to   estimate   x.    and   c.    in   conjunction 
with   the  dimensionless   terms   defined    in   Eq.    11.      If  the  values   of 
E,     lie   outside   the   range   shown   on    Fig.    6  then    Eqs.    7   and   8  should 
be   used   to   calculate   x.    and   c,  ,    respectively.      The  use   of  Fig.    6 
in  design   computations    is   described   later. 


2.3.3         Effects  of   Neglecting   Side-Wall    Reflection   of  Material: 

As    stated   previously,    the   critical    co-ordinates   pr-edicted  by   Eqs.    7 
and   8  are   subject   to    some  limitations   due   to   the   assumption   of  n=0. 
An  empirical    procedure,   described    in   previous    publications    (7,   8), 
was    utilized   to   determine   the   region   of  mixing   zone  within   which   the 
predicted   pollutant   concentrations   are   not    likely   to   be   affected   by 
the   reflection   of  the  material    from   the   confining   side-walls   of  a 
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FIGURE  6.   NON-DIMENSIONAL  CO-ORDINATES  OF  CRITICAL  POINT 


river  channel.  Using  those  results,  a  graph  of  log  9  versus  p 

m 

was   prepared    (Fig.    7),  where   p^  denotes  the  maximum  dimensionless 

value  of  p   at  which   the   total    number  of  images  was   zero.      The  values 

of  Cp  at   all    points  which   lie   under   the  curve  are  not   likely  to   be 
affected  by   the   side-wall    reflections. 

By  making   use  of  Fig.    7  and   the   values   of  p.  ,   E.    and   G.  ,    it   is 
possible   to   examine  whether  the   concentration   at   the   critical    point 
is   affected   by  the   side-wall    reflections   or   not,   as   follows:      From 
Fig.    7,  we  can   obtain   the   value  of    v    corresponding  to   p     =  p.; 
the  value  of    9    ,    can   be   calculated   from  Eq.   27.      If    cp    .    is 
less  than    9  ,   then  the  effects  of  side-wall    reflections   at  the 
critical    point   are  negligible,  whereas    cp      greater  than    9 
indicates   that   there  may   be   some   approximations    in  the  estimation  of 
c.    and   c    ..      In   such   cases.   Fig.    4  can   be   used   to   predict   the 
longitudinal    concentration   distribution   along  the  known   lateral 
boundary  of   LUZ;   the   value  of  c,    is   then   obtained    from  this   curv( 


/e, 


2.3.4    Maximum  Longitudinal  Spread 

In  the  case  of  a  bank  outfall,  the  maximum  longitudinal  spread  of  a 
LUZ,  X   ,  associated  with  an  effluent  concentration,  c  ,  may  be 
obtained  by  plotting  c(x,o)  vs.  x  and  drawing  a  horizontal  line 
through  c  to  intersect  this  curve  at  A,,  as  shown  on  Fig.  8. 
The  concentration  distribution  along  a  lateral  boundary  of  LUZ, 
c(x,  q.  ),  is  shown  by  the  curve  2  in  Fig.  8;  from  this  curve,  a 
critical  concentration,  c.  ,  is  obtained.  (Note:  The  critical 
concentration  can  also  be  determined  directly  from  Eqs.  7  and  8,  as 
out  1 ined  earl ier) . 

When  the  effluent  concentration  is  equal  to  C  .  (determined  from 

Eq.  6),  the  concentration  distributions  along  the  shoreline  and 

along  the  lateral  boundary  will  be  those  indicated  by  the  curves  3 

and  4,respecti vel y  (see  Fig.  8).  Then,  the  length  x.^^., 

associated  with  c  .,  is  obtained  by  drawing  a  horizontal  line 

through  c^  to  intersect  the  curve  3.   It  is  interesting  to  note 

that  X   .  can  also  be  obtained  by  drawing  a  horizontal  line 
sceA  -^ 
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FIGURE  7.   SIDE-WALL  REFLECTION   EFFECTS 


N 

D 

J 

X 

0 

^ 

r— 1 

Q 

^ 

< 

(T! 

W 

14-1 

a 

-P 

cu 

a 

CO 

0 

J 

e 

< 

0 

2: 

V-l 

*-H 

•4-1 

Q 

0) 

EH 

u 

i-i 

c 

■n 

03 

2 

-U 

0 

CO 

J 

•H 

T3 

0 

E 

&-, 

(T3 

0) 

Q 

i-l 

0 

4J 

X 

to 

Eh 

c 

U 

s 

2 

0 

Q 

< 

M 

< 

n 

CM  d) 

cu 

-St   o 

<: 

en 

a; 

X 

0 

CO 

LU 


{"^h'x)  D    '  {o'x)o 


-    24    - 


through  c.  to  intersect  curve  1  at  B,  and  then  drawing  a 

vertical  line  BiB-.  The  latter  method  of  determining  x   , 
12  ^  sceA 

will  circumvent  the  need  to  plot  curves  3  and  4  when  the  lateral 
boundary  and  all  other  design  parameters  remain  unchanged. 

In  some  cases,  plots  of  c(x,o)  vs.  x  may  result  in  curves  which  are 
almost  horizontal  at  large  values  of  X,  thus  making  it  somewhat 
difficult  to  determine  the  intersection  points  A,  and  B, 
accurately  by  the  graphical  procedure.  The  use  of  MIXAPPLN  program, 
described  in  Section  2.6,  will  provide  more  accurate  values  of  the 
longitudinal  spreads. 

2.4  Effects  of  Backgrouno  Concentrations 

The  foregoing  formulations  have  been  based  on  the  assumption  that 
the  concentration  of  the  nonconservati ve  pollutant  of  concern  in  the 
nver  just  upstream  of  the  outfall  is  zero.  When  the  background 
concentrations  are  not  negligible,  it  is  possible  to  include  their 
effects  by  assuming  that  the  background  concentration  distribution 
in  the  cross  section  of  the  river  is  uniform  throughout  the  study 
stretch.  This  assumption  is  valid  in  cases  where  the  effect  of  any 
upstream  discharge  of  pollutants  results  in  a  uniform  cross- 
sectional  distribution  just  upstream  of  the  effluent  discharge  under 
consideration.  The  decay  rate,  k,  ,  of  background  pollutants  may 
differ  from  that  associated  with  the  effluent  discharge.  The 
concentration,  c.  ■ ,  at  each  transect,  i,  due  to  an  upstream 
background  value,  c.  ,  is  calculated  from 


c 


bi  =  S  exp(-k^^x./u^)  (12) 


in  which   u.    =  average  value  of  velocity   from  the  outfall    to   the 
transect    i  ;   and    k   -j.   =  temperature-corrected   value   of   k,  .      The 
value,   c.     ,    is   then  added   to   each  concentration   value,   c(x,   q),   to 
give   the   total    concentration    resulting    from   the   upstream   and 
effluent  discharges   of  the   pollutant. 
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When   there  are   a   series   of  outfalls  discharging   into   the   same   river, 

the  effect  of   an   upstream  discharge  just   above   the   next   downstream 

outfall    would   depend   on   the   distance   between   the   two  outfalls,   x    , 
'^  w 

as  well    as   on   the  characteristics  of  the  waste   constituent   and   the 
hydraulic   characteristics   of  the   river   channel.     When   the   oufalls 
are   sufficiently    far   apart,    the   cross-sectional    concentration 
distribution   would   likely  approa:h   uniformity.      In   such   cases,   the 
one-dimensional    assumption   could   be   used   to   estimate   concentrations 
in   the   river  just   above  the  downstream  outfall    locations    so  that 
these  values   can   be   used   as   background   levels    in  modelling   the 
latter   discharge. 

2.5  Scale-up  of  Model    Parameters 

When   the  model    is   applied   to   predict   pollutant   concentrations    in 
mixing  zones   for   specified  managtiment  options   such   as   treatment 
plant   expansion,    flow  augmentation,   seasonal    variabilities    in 
streamflow  and   temperature,    etc.,    it    is   necessary   to   scale   up  the 
model    parameters   to   account   for   the   variabilities   associated  with   a 
defined  option.      The   relationships   commonly   used   to   scale   up  the 
parameters   affected   by  changes    in   design   temperature  and   design 
streamflow   are  presented   below. 

The  temperature  dependence  of  the  decay   rates    is  expressed   by  the 
modified  van't   Hoff -Arrhenius    relationship: 


e  "^2   -  ^1 

K-|.„-K-|-i  \  '■■^  1 


in   which   k,,    and   k-|.„   =  decay   rate   coefficients   at   temperatures 
t'  C  and   T^C,    respectively;   and      e  =  temperature   correction 
factor.      The  values  of   e  for   TRC   and   aimionia  may   be   taken    as   1.03 
and    1.106,    respectively   (3,    4).      The  decay   rates  may   also  depend   on 
turbulence   and   other   hydraulic   characteristics   of   stream  channel, 
and   effluent   and   stream  water   quality;    however,   only  the   effect   of 
the   primary   variable   temperatu'-e ,    is   considered    in   this   study. 
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A  change  in  river  discharge  results  in  corresponding  changes  in 
depth,  H,  velocity,  U,  and  top  width  of  water  surface,  B.   In  water 
quality  studies,  the  changes  in  the  hydraulic  parameters  are 
generally  expressed  by  the  Leopold-Maddock  equations  (3) 


B^  =  B^  (Q,/Q2)''^'  (14) 

H^  =  H^  (Qt/Q2)^^'  (15) 

U^  =  U^  {Q^/Q^f^^  (16) 


in  which  B  ,  H  and  U  represent  in  order,  the  mean  values  of 
m   m     m 

width,  depth  and  velocity  at  a  flow  Q  (m  =  1 ,  2);  and  bex,  hex, 

and  uex  are  exponents,  such  that  their  sum  is  equal  to  unity.  The 

values  of  the  exponents  are  estimated  from  the  channel  geometry, 
streamflow  and  time  of  travel  data. 


The  effect  of  a  change  in  streamflow  rate  on  the  dimensionl ess 
coefficient,  S  ,  is  not  known;  in  this  study,  the  coefficient  is 
assumed  to  be  independent  of  discharge  in  a  given  river  stretch. 

2.6  Computer  Program  MIXAPPLN 

A  computer  program  MIXAPPLN  (Mixing  zone  -  APPLicatioN  model), 
written  in  FORTRAN  language,  has  been  developed  for  performing 
design  option  computations.  This  program  has  been  documented  in 
Water  Resources  Paper  #14  (Ref.  8).  The  program  input  includes 
essentially  the  same  data  as  the  MIXCALBN  model,  as  well  as  data  on 
options  normally  consider^ed  in  design  runs  (viz..  numbers  and  values 
of  upstream  flow  rates,  effluent  flow  rates  and  temperatures,  as 
well  as  pH  values  in  the  case  of  ammonia). 

The  program  MIXAPPLN  is  set  up  to  handle  eight  upstream  flow  rates, 
six  effluent  flow  rates,  and  six  temperature  values,  as  well  as  four 
pH  values  in  the  optional  case  of  prediction  of  toxic  airmonia 
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levels.  For  each  run,  the  program  output  consists  of  a  table  of 

transverse  concentration  distributions  at  each  transect  for  p.  =  0 

to  0.5  and  a  second  table  which  includes  c,  .  x.  ,  c  .  and 

X   „  for  p,  =  0.2,  0.3  and  0.4.  The  maximum  longitudinal 
sceA    "^L 

spread,  x   ,  corresponding  to  the  input  effluent  concentration, 

c  ,  is  also  printed  out.   Some  of  the  results  of  each  run  are 

e  ^ 

stored    in   arrays   and   printed  out   as  a    surmary   table. 
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III.    DATA  COLLECTION   AND   ANALYSIS    PROCEDURES 

3.1      Field   Survey   Procedures 

The    field    surveys   are  designed   to   collect   data   on   outfall    discharge 
and   background  water   quality   characteristics,    and   on   the   transverse 
distribution   of   various   parameters   of   interest   at    four  to    six 
transects.      The  details  of    field   surveys   are   usually   dependent   on 
site-specific   conditions.     A  general    field   study  procedure   is 
described   below. 

The   location   of  transects  can   be   based   on   preliminary   in   situ 
measurements  of  a   conservative   parameter    (e.g.    conductivity)    at 
selected   access   points   to   establish   the   approximate   longitudinal 
boundary  of  the  mixing   zone.      The   selected   locations   of   transects 
must   be  marked   on   a  map   so  that   they  can   be   readily   identified    in 
the   field  during    further   surveys.      Cross-sectional    depths  must   be 
measured   at   a  minimum  of   15  points   at   known   lateral    distances 
measured   from  a    reference   bank    (usually   the  outfall    bank).      Velocity 
measurements  must   also  be   taken   at    least   at   two  transects   by 
following    standard    streamflow  gauging   procedures;    however 
measurements   at   all    transects   are   desirable.      (Note:     At   transects 
where  velocities   are   not  measured,    the  Manning's  equation   can   be 
used   to   simulate   the   velocity   profiles   using  measured   depth   profiles 
as  outl ined   1 ater) . 

Water   samples   must   be   collected   at   the   upstream  boundary  and   at   the 
effluent  outfall,    and   at   each   point   where   cross-sectional    depths   and 
velocities   are  measured.      However,   the   samples   can   be   collected   at 
selected   points   at   each   transect    (viz.,    less   points  outside  effluent 
plume,   alternate   points,   etc.)    to   reduce   sample   analysis   costs,    in 
which   case   the   concentrations   at    other   points   are   obtained   by 
interpolation   during  data   analysis.      In   order   to   account    for   the 
possible   effects  of  fluctuations    in   effluent  water   quality   and 
discharge   on   the    instream  concentrations,   the   sampling    is   carried 
out   either  by   following   the    same  plug  of  water   beginning   at    the 
outfall    and   proceeding   to   successive  downstream  transects,   or   'iuring 
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a    round-the-clock   intensive  survey  when   samples   are  collected   at 
each  point   at   3   to  4   hr.    intervals.      Obviously,    the   selection   of  a 
sampling  methodology  would   depend   on   the  man-power,   time  and   other 
resource   constraints,    as  well    as   the   objectives  of   the   study. 

In   situ  measurements  of  temperature,   pH  and   conductivity  must   be 
taken   along  with   collection  of   samples  which   are  to   be   analysed    in 
tlie   laboratory  for   nonconservati  ve  pollutants   of  concern.      In   some 
cases,    it  may   be   desirable   to    inject  a   solution   of  dye   continuously 
to   gather   data   on   the   transverse   distribution   characteristics  of  the 
river.      This    is   particularly   useful    to   simulate  effluent  discharge 
from  proposed   outfall    locations,   and   in   cases  where   relocation   of  an 
existing  outfall    is   being   considered.      The  dye    injection  must   be 
maintained    for   2  to   3  hours   or  more  to   establish   steady  state 
conditions.      The   cross-sectional    distributions  of  dye   at   selected 
transects   can   be   obtained   directly  through  fluorometric   tracing; 
however,    in   shallow   rivers  where   the   passage  of  a   boat    is   difficult, 
it  will    be   necessary  to   collect  water   samples   at   known  cross- 
sectional    points   for  dye   concentration   analysis. 

Generally,    two  surveys   must   be   carried   out   under  different   instream 
hydraulic   conditions   so   that   the  data  of  one   survey   can   be   used   to 
calibrate   the  model,   which    is   then   verified   using   the   data   of  the 
second   surv;?y.      If  man-power   and   other   resource   constraints   dra 
limited,   then   the  data   of  one   survey  may   be   used   to   validate   the 
model.      The   limitations  of   the   letter   approach  must   be  given  due 
consideration    in   using   the   predictions   of  the  model    for   planning   and 
management   purposes. 

3.2  Data   Analysis   Procedure 

The  data   collected   during   one   of  the    field    surveys   are   utilized   to 
determine   the   parameters    required   for  modelling.      A  computer   program 
MIXANDAT,   written    in    FORTRAN   language,    is   utilized   to   perform   the 
following   computations   using   the    survey   data   as    input: 

1.        Average  depth   and   velocity   at   each   transect. 
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2.  Simulation   of   velocity   distributions   at   cross-sections 
where  velocities   are   not  measur^ed. 

3.  Transverse   distribution   of  streamflow  at  each  transect, 

4.  Mass   flux  values   of   conservative   and   nonconservati  ve 
materials  at   each   transect. 

The   program  output   also    includes   shape-velocity   factors   and  variance 
of  cross-sectional    distributions   of  conservative  materials   (required 
to  evaluate   the  dispersion   characteristics);   the   pertinent 
computational    details  are   presented   in   Reference  3.      A  simplified 
method  of  evaluating   the  dimensionl ess   coefficient,      s    ,    is 
presented   later.      The  computations    involving  mass   flux  and    variance 
utilize   net   concentration  values  of  a  material    determined   from  the 
measured   values   minus   background   value  of  that  material    in    river 
water.      If  a   net   concentration    is   negative,    then    it    is    set  equal    to 
zero. 

A  brief  description   of  the   computer   program  MIXANDAT   is   presented   in 
Section   3.4.      All    calculations    in   the  MIXANDAT  program   are   based   on 
the   trapezoidal    approximation   of  the   integrand.      Details  of  the 
various   computations   are   presented   below   (with   the  exception   of  the 
variance  computations).      For   simplicity,    the   subscript,    i, 
associated  with   the   transects  will    be  dropped   from  the 
relationships,   except  where   it   is   necessary  to   retain    it   for   clarity. 

3.2.1  Mean   Depth   and   Velocity: 

At  a  given   transect,    i,   the  cross-sectional    average  values   of  depth, 
h,   and   velocity,    u,    are   computed    from: 

R   =  A/b  (17a) 

Û   =   Q/A  (17b) 

in   which  A   =   cross-sectional    area   of   flow;    and   b   =   top  width  of 
water   surface   at   the   transect,    i.      The   cross-sectional    area    is 
computed    from: 
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m 

A  =     I   0.5(h^   +  ^.i)(^kn   -^k)  (18) 

k  =  l 

in   which  m  =  number   of   vertical    strips    in   the   cross-section;   and 
h,     =  local    depth   at   a   lateral    distance,  y,  ,    from  the   reference 
bank   of   the    river.      The   total    discharge    in   the   river,    Q,    is   given  by 
the   sum  of  upstream  and   effluent   flow  rates. 


The  average  values  of  width,   depth   and   velocity  within   a  given 
channel    stretcl 
rel ationships : 


channel    stretch  of  length,   x    ,   can   be   computed   from  the   following 


m-1 


1  =0 


(b^    -b.^^)(x.,^-x.) 


B     =   1/2 


^.1 


(19a) 


m-1 

z 

i=o 


(^    ^^.l)(^    ^°i.l)(^.l  -i) 


H      =      1/4 


X  . 


i+1 


(19b) 


U      = 


BH 


(19c) 


where   Q,    B,    H  and    U  denote   n  ver   discharge,   average  channel    width, 
mean   depth   and  average  velocity,    respectively. 

3.2.2         Transverse   Distribution   of   Flow: 

The   cross-sectional    distribution   of   streamflow   at   a   given   transect 
is  calculated    from 


\  =o-25(\  ^^.i)^^  ^^.i^^^k.i  -^k) 


:20} 
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J 

QJ      =        ^       q,.  (21) 


in   which   q,     -  discharge   through  an   elemental    strip,    k ;    Q' .    = 
partial    cumu'ative  discharge  between   the   reference  bank   and  j 
strip;   and   u,    =  depth-averaged    local    velocity   at   a    lateral 
distance,  y,  ,    from  the   reference   bank.     Measured   velocity   values 
are   used    in   Eq.    20  whenever   such  data   are  available.     At   transects 
where  depth   profiles   are  available   and   velocity  measurements   are   not 
taken,    the  velocity  distributions   are   synthesized  by  making   use  of  a 
resistance   equation    (either   Manning's   or   Chezy's   equation)    following 
Pal  ley   and   Sayre    (11),    as  outlined   below. 


In   a   shallow,    wide   river   the   local    velocity   u,   at   a   point   where   the 
depth   IS  h,    can    be   expressed  by   the  Manning's   equation 


1      2/3     1/2 
u^.h         s^  (22) 

in  which   n'    =  the   Manning's   roughness  coefficient;   and   s     =  energy 
gradient    (approximated  by   the   slope  of   the  channel    bed   surface). 
Similarly,    the   cross-sectional    mean   velocity   is   given   by: 


^     2/3      1/2 
'-n'^  'e  (23; 

From  Eqs.    22   and    23  we  get: 


;u/u")  =  (h/h)2/^  (24; 


This  expression   can   be   generalized   to   account   for  possible  variation 
of  bed    roughness    in   a   given    stretch,   as    follows: 


-   ^2 
;u/u)    =   f^    (h/h)    \  (25) 
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in  which   f,    and    f^  are   empirical    coefficients.     According   to 

Manning's   equaion   f,    =1.0  and    fp   =  2/3;   when   Chezy's  equation 

is   used    f,.    =  1.0  and   fp   =  1/2.      The   values   of   f,    and   f„  can 

be  determined  by    fitting   the  velocity  distribution  data   collected   at 

a   few  transects   to   Eq.    :15.      In   the   program  MIXANDAT,   the  default 

values    for   f,    and    f^  art;   1.0  and   0.67,    respectively. 

Eq.    25    is   used   to   simulate   the  velocity  profile   at   each   transect 
where   velocities   are   not  measured    (provided  depth   profiles   are 
available).      The   simulated   or  measured  velocities   are   used   to 
compute   a  value  of  total    discharge  Q' ,   at   a   given   transect   using 
Eq.    21.      If  Q'    differs   from  the  known   total    discharge,    Q,    the 
velocities   are  multiplied   by  a  correction    factor,    (Q/Q'),   to   obtain 
a   new   set  of   simulated  values   at   each   transect,   and   the   latter 
values   are  used    in    further  computations. 

3.2.3         Mass   Flux   and   Relative  Concentrations: 

The  mass   flux,    f,  ,    of  each  material    (conservative  or 
nonconservati ve)    passing   through  a   vertical    strip,    k,    at   a  given 
transect,    i,   is  determined   from   the   relation 

^k  =  n/8)(y,,i    -y,)(h^  ^h,^^)(u,  -u,^^)(c,  .c,^^),  (26) 

in  which  c,  represents  the  measured  concentration  at  y,  ,  minus 
background  value.  The  total  net  flux,  F.,  of  each  material  at  the 
transect,  i,  is  then  determined  from 

m 

F,  =  ^   f,  (27) 

k  =  l 

The   total    gross  flux  of  a  conservative  material    at   each   transect 

(equal    to   total  net   flux   plus   background   flux)   must   be   equal    to   the 

sum  of  effluent  and    upstream   flux   values.      However,    some  deviations 

generally   occur  in   pr.^ctice  due    to  various    reasons   which   are 

site-specific.      The   cross-sectional    average   concentration,   c   .     , 

a  t  r 
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of  each  material  is  computed  from  c,.   =  F./Q.  Nondimensional 

at  r         1 

(i.e.,    relative)    concentration   values  of   conservative  materials   are 

then   calculated    from  c       =  c,  /c  .     ;   these   values  will    be   used 

ra         K     atr 

for   comparison  with   the  predictions   of   the   stream  tube  model.      The 

use   of  relatii/e  concentration   values    removes,   to   some  extant,   the 

effects  of   deviations    in   the  mass  flux  values  of   conservative 

materials  at   various   transects,   thus   mîking   c       values   better 

r  a 

suited   for  comparative   purposes  during  model    validation   studies. 


The   total    gross   mass    flux   values   of  a   nonconservati ve  material    at 
various   transects   are   used   to  determine   the   first-order  decay   rate 
of  the  material    by  preparing   a   semi-log   plot   of  mass    flux   versus 
time  of  travel.      (Note:    Sometimes,   distance   is   used    in   plotting 
instead  of  time  of  travel  ,    in   which   case  mean   velocities   are  also 
required   to   compute   the  decay   rates.) 

The   units   of  the  mass    flux  values   determined   from  Eq.    26   and   27   are 

dependent   on   the   units   of  various   quantities   appea^^ing   on   the   right 

hand    side   of   Eq.    26.      Thus,    if   length    is   in   meters,   time   in    seconds 

and   concentration    in   ug/L,   then   the   units  of  mass   flux  will    be   1000 

3 
ug  per   second   or  mg/s    (since   1    m     =  1000  L);    if  the  concentrations 

are    in   mg/L,    then   the  mass   flux   units   are   gm/s. 


3.3 


Simplified  Method  for  Estimation  of  s 


A  simple  method  for  the  determination  of  S  can  be  developed  from  the 

following  equation  for  the  peak  concentration  of  a  conservative 

tracer,  obtained  from  Eq.  2  with  n  =  0  and  k,  =  0: 

d 


'^''  °)  ^  w  =  s(-"-^  ) 


-1/2 


(28) 


in  which  c    =  maximum  concentration  at  a  distance  x  below  the 
max 

outfall  and  c  =  c  Q  /Q.   Substituting  cp  =  Sx/B  and 
simpl 1 fy ing,  we  get  : 


^'a/c^ax)  -   (  "^  )(x/b) 


(29) 
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Thus,  a  plot  of  (c  /c   )^  versus  (x/b)  would  result  in  a 

a  max 

straight  line  of  slope,  s   =   {  ^s,    );  then  B  =  s/t^  . 


This  method  is  much  simpler  than  the  variance  method,  requiring 
channel  widths  and  peak  concentrations  at  various  transects  below 
the  outfall,  and  values  of  c^,  Q  and  Q.   It  shou:d  however,  be 
noted  that  the  method  gives  an  approximate  value  suitable  as  a  first 
estimate  for  modelling  studies. 

3.4     Computer  Program  for  Data  Analysis 

As  stated  previously  the  computer  program  MIXANDAT  (Mixing  zone  - 
ANalysis  of  DATa)  is  used  to  analyse  the  data  collected  during  field 
surveys.  The  program  is  set  up  to  accept  data  at  up  to  eight 
transects  and  up  to  50  values  of  lateral  distances,  depths, 
velocities  and  concentrations  at  each  transect,  as  well  as  up  to 
five  water  quality  constituents  at  each  transect.  A  table  of 
widths,  depths,  velocities,  concentrations  and  mass  flux  values  is 
printed  out  at  each  transect  for  each  of  the  constituents. 

The  variance  values  at  each  transect,  computed  from  different 
methods  by  considering  the  distribution  of  each  water  quality 
constituent,  are   also  tabulated.  The  program  is  documented  in  the 
Water  Resources  Paper  #14  (Ref.  8). 
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IV.  APPLICATION  OF  PREDICTION  TECHNIQUES 

4 . 1     General  Procedure 

The  general  procedure  involved  in  the  application  of  the  assessment 
technique  for  mixing  zones  can  be  sumnarized  as  follows: 

1.  Preliminary  study  to  gather  information  on:  streamflow, 
temperature,  pH,  average  channel  width,  depth  and  velocity,  type 
of  outfall  and  its  location  in  river  channel,  effluent  flow  rate 
and  concentration  of  pollutant,  tracer  distributions  (eg. 
conductivity)  at  easily  accessible  locations  in  the  river 
channel ,  etc. 

2.  Preliminary  predictions  of  cross-sectional  concentration 

distributions  at  selected  downstream  locations  using  estimated 

val ues  of  s  and  k , . 
d 

3.  Selection   of  the   locations   of  transects    for   detailed   data 
collection,    based   on   the  prel  im  i  nai'y    field  measurements   and 
predictions  . 

4.  Collection   of   field   data   as   outlined    in   Chapter    III. 

5.  Analysis  of  data   and   estimation   of   parameters    (see  Chapter    III 
for  detai 1 s) . 

6.  Model    validation   studies   to   determine   the  model    parameters     £  and 
k,    (see  case   studies   in   Water   Resources   Paper   #14). 

7.  Evaluation   of   the   impact   of   various  management   options   on   the 
mixing    zone   boundary  and   allowable   effluent   concentration   of 
pollutant  by  making   use  of   the   step-by-step  design   procedure   for 
critical    point   method   outlined    in    Section    4.2,    or   by  using   the 
computer   program  MIXAPPLN. 

8.  Specifying   acceptable  mixing   zone   boundary  and    associated 
permissible   effluent   r.oocent'-'at  i  on  , 
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4.2  Design   Procedure   for  the  Critical    Point   Method 

A   14-step  design   procedure    Is   developed   for   the   application   of  the 

critical    point  method.      This   procedure   is   based   on   the   use   of  the 

graphical    solutions    (Figures   4  and  6)   presented   in  Chapter    II.      It 

consists   of  three  major   parts:    Steps    1  to   4  involve  the   selection   of 

model    parameters   for  a   design   case;      steps   5  to   11   are   related   to 

the   computation   of   x,  ,   c,    and   c   „;   and   the  determination   of 
^  L       L  eA' 

longitudinal  boundary  of  LUZ  is  outlined  in  steps  12  to  14.  An 
example  which  illustrates  this  design  procedure  is  presented  in 
Section   4.3.      The   step-by-step   design   procedure    is   as   follows: 


1.      Select   parameters    for  the  design   condition,   viz.,    Q,    Q   ,    c 
temperature   and   pH. 


2.      Determine  the   channel    hydraulic   parameters   B  and   U,   and   the 
coefficient,   8    . 


3.      Obtain   the   pollutant  decay   rate,   k  . ,   and   apply  temperature 
correction,    if   necessary. 


4.      Determine   allowable   instream  concentration,   c   ,    In   accordance 
with   the   Provincial    Water   Quality  Objective   (PWQO)    for   the 
pollutant   under   consideration    (Ref.    10). 

::U-f  •  5.      Calculate   the   cross-sectional    average  concentration    from 


.<. 


c     =  c   Q  /Q. 
a         e^e   ^ 


6.   Select  the  lateral  boundary  of  LUZ,  q.  ,  and  calculate  p.  = 

q./Q  (from  Eq.  11).  Obtain  the  value  of  ?  corresponding  to 

p  =  p,  from  Fig.  7. 
•^m   ^L       ^ 


7.      Calculate   E,    from   E^^   =  k^B/8U, 


8.      From   Fig.    6,   determine   the   values   of  G,    and   c     corresponding 
to    E^   and   p.    obtained    in    steps   6  and   7.      If   the   values   of 
E,    and   p.    àrti   outside   the   range   shown    in    Fig.    5,    use   Eqs .    7 
and   8  to   compute   x.    and   c,  ,    and   proceed    to   step    10. 
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9. 


Calculate   the   value   of     cp   .    =  Gi /E,     from   Eq.    11).      If   this 

value    is   less   than   or  equal    to   the  value  of     cp  obtained    in   step   5, 

then   calculate   the  co-ordinates   of  the   critical    point,   x,    and 


c,  ,  from  X.  =  G,  U/k  ,  ,  and  c. 


^r^a 


Proceed  to  step  11, 


10.  If  cp  ,  is  greater  than  cp  ,  use  Fig.  4  to  calculate  concentrations 
at  several  values  of  x  for  the  selected  p.  and  prepare  a  plot 
of  c(x,  q.  )  vs.  x;  determine  c.  from  this  plot. 


11.  If  c,  is  not  equal  to  c  ,  then  calculate  c  .  from  Eq.  6. 

L  s  eA      ^        ■  '  ■ 

This  value  represents  the  allowable  effluent  concentration  based 

on  the  LUZ  concept.   If  c,  =  c  ,  then  c  =  c  ..   Compute 

L    s'      e    eA     '^ 


percent  treatment  =  100  (c 
than  c  . 


c^^)/c^,  when  Cg^  is 


less 


12.  Calculate  the  pollutant  concentrations,  c(x,o),  along  the 

discharge  shoreline  for  several  values  of  x,  by  making  use  of 
Fig.  4  (as  per  the  procedure  described  in  subsection  2.1.3)  and 
plot  c(x,o)  vs.  x. 


13.  Draw  a  horizontal  line  through  the  instream  pollutant 
concentration  objective,  c  ,  to  intersect  the  shoreline 
concentration  profile  at  a  point  A^  and  draw  a  vertical 
A.A„  (see  Fig. 8).  Then,  the  spread  resulting  from  the 


effluent  concentration  c  is 

e 


X    =  0A„. 
see     2 


1  ne 


14.  In  order  to  determine  the  spread  x   .  due  to  the  allowable 

sceA 

effluent  concentration,  c  .,  draw  a  horizontal  line  through 
c.  to  intersect  the  shoreline  concentration  profile  at  a  point 
Bl  and  draw  the  vertical  line  B.Bp  as  shown  in  Fig.  8. 
Then,  x^^^^  =  OB^. 

4 . 3  Total  Ammonia  in  the  Mixing  Zone  of  Grand  River  below  Water! oo 

In  order  to  illustrate  the  computational  aspects  of  the  methodology, 
examples  based  on  the  results  of  studies  in  the  Grand  River  below 
the  Waterloo  pollution  control  plant  are   presented  in  this  section. 
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These   studies  are  described    in   other   publications    (12,    17)   and 
should   be    referred   to   for  more   details. 

4.3.1         Illustration   of   Step-by-Step   Design    Procedure: 

An   example   based   on  a   survey   carried  out  during  June   1975    in   the 
Grand   River   below  the   Waterloo   pollution   control    plant   outfall    has 
been    selected   to   illustrate   the   14-step   design   procedure.      Using   the 
chloride   ion   concentration   distribution   data,  the   value   of  B     was 
estimated   to   be   5.66   x   10'      (from  the   variance  method);   however, 
from  model    validation    studies,    the   coefficient  was   found   to   be   a 

function   of   longitudinal    distance   below  the  outfall;   its   value 

-3  -3 

ranged   from  0.8  x   10       to   2.5  x   10     .      For   the   purpose   of 

-3 

illustrating   the   step-by-step   design   procedure,      B    =   1.0  x   10 

has   been    selected   to   be   the   representative  value.      The  effluent 
concentration   of  total    arrmonia   used    in   the  computations    is   12.5  mg/L 
(as   N)    and    is  typical    of  a   conventional    activated   sludge  plant.      The 
decay   rate  of  total    aixionia   v;as   obtained    from  the   results   of 
wastewater   assimilation    studies   carried  out    in    tiiat    stretch  of  the 
ri  ver . 

Design    Problem:    Total    AiTinonia    in   the   Mixing   Zone  of  the   Grand   River 
below  Waterl oo. 


1.  Q  =   12.54  m^/s;     Q^   =  0.39  m^/s    ; 

c      =   12.5  mg/L   (as    N);      Temperature   =   20°C;      pH   =  8.3. 


57.32  m;        U   =  0.39  m/s;        S    =   1.0  x   10 
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3.       k     =  1.71  X   10"^/s  =   1.48/day  (to  base  e)   at   20°C. 


4.        From   Ref.    10,    the   PWQO   for   un-ionized  ammonia    is   0.02  mg/L 
(as    N);      in   terms  of   total    armionia,     c.    =   0.272  mg/L   (as 
N)    fi-om   Eq.5. 


5.        c     =    (12.5)    (0.39)/(12.54);      =  0.389  mg/L 
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6.        q^   =  5.015  m-^/s    ;        p^   =   (5.016/12.54)    =   0.4; 
from  Fia.    7,     cp    =   0.065. 


7.        E^   =   (1.71   X    10'^)    (57.32)/(10"^)    (0.39);      =     2.51 


From  Fig.    6,      G,     =  0.15   and   c     =   1.0. 


9.        cp      =   (0.15/2.51)    =   0.06;   this    is   less   than     cp   =  0.065, 
x^   =   (0.15)    (0.39)/(1.71   X   10'^);      =  3421   m. 
c      =    (1.0)    (0.389)      =  0.389  mg/L. 


10.      (This   step    is   not   applicable] 


11.      c  ^  =   (12.5)    (0.272)/(0.389);      =  8.74  mg/L   (as   N) 
%  treatment   =    (  100) ( 12. 5-8. 74)/( 12.5) ;      =   30.1% 


12.      Table   2  shows   the   computation   of  c(x,o)   using   Fig.    4   (see 
subsection   2.1.3).      A  plot   of  c(x  o)    versus  x    is   presented 
in    Fig.    9. 

TABLE    2.      SHORELINE   CONCENTRATION   DISTRIBUTION   -   DESIGN   PROBLEM 


X 

"P 

Cp 

Cc(x,o) 

c(x,o) 

m 

mg/L 

mg/L 

200 

0.0035 

9.5 

3.70 

3.67 

500 

0.0087 

5.0 

2.33 

2.28 

1000 

0.017 

4.2 

1.63 

1.56 

2000 

0.035 

3.0 

1.17 

1.07 

4000 

0.070 

2.0 

0.78 

0.65 

5000 

0.10 

1.7 

0.66 

0.51 

10000 

0.17 

1.3 

0.51 

0.33 

15000 

0.26 

1.2 

0.47 

0.24 

20000 

0.35 

1.1 

0.43 

0.18 

13.  From  Fig.    9,   the   longitudinal    boundary  of   LUZ   resulting 
from  the   effluent   concentration   c      is:   x  =   12800  m. 

14.  The   spread   due   to   allowable   effluent   concentration   c    . 

is:      X        „    =  8600  m,    (obtained   by  drawing   a   horizontal 
sceA  '  -^ 

line   through  c^^   =  0.389  mg/L). 
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4.3.2.       Management   Options: 

The   impact   of  several    management   options   on   the  mixing   zone  boundary 
and   allowable   effluent   concentrations  of   total    ammonia  were 
predicted   using   the   computer   program  MIXAPPLN.      The  options 
considered    in   this   study   are   as   follows: 

(i)         Two  effluent   flow  rates,    0.307   and   0.60  m  /s, 

representing   the  present   and  a   future   hydraulic  capacity  of 
the  plant. 


(ii)  Two  river  discharge  values    (just   upstream  of  the   outfall), 

3 
4.24  and   8.49  m  /sec,   attainable  without   and   with   flow 

augmentation   from  existing   reservoirs. 


(iii)  Three   temperature   values,    15    ,   20     and    25   C, 

reflecting   the    range  of   seasonal    variability  during   the 
open  water  season    (May-October). 


(iv)         Two   pH  values,    7.8  and   8.3  units,   which   are   also   based   on 
observed  seasonal    variability  of  this  parameter   in   the 
Grand  River. 

There  are   24  combinations   of  the   parameters   associated   with  the 
above  options,    requiring   24   separate  runs  of   the  model.      The  options 
listed   above  are   extracted    from  a  technical    report   prepared   by  Post 
and   Gowda    (12),   which   deals  with   effluent  mixing   zone   studies 
undertaken   as   part   of  the   Grand   River   Basin   Water  Management   Study 
(GRBWMS). 

The   channel    hydraulic   parameters    (viz.    river  discharge,   width, 
average  depth   and  mean  velocity)    obtained   from   the  June   1975   field 
test,  were   used   as    reference   values.     The  width,  depth   and   velocity 
parameters  were  scaled-up   for   the  design   streamflow  conditions   for 
each   option    listed   above  using   the   Leopol  d-Maddock   equations    (Eqs. 
38-40)    in   which   the  values   of   the   exponents  were:      bex   =   0.05,    hex    = 
0.50   and    uex   =   0.45.      The   background   concentration   of  total    anmonia 
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in   the   river    (just   upstream  of  the   Waterloo  WPCP  outfall)   was 
assumed   to   be  zero.      The   Kitchener  WPCP  outfall,    located   20.487 
meters   below  the   Waterloo  outfall,   was   considered   to   be   the   limiting 
downstream   boundary. 

The  concentration   of  total    amnonia    in   the    Waterloo  WPCP  effluent    is 

taken   to   be    12.5  mg/L    (as   N),    as    stated   before.      The    instream  decay 

-5 
rate  of  total    anmonia  was   taken   as    1.71   x   10     /sec   or   1.48/day   (to 

base   e)   at   20  C,    based   on   previous  waste  assimilation   survey 

results;   the  decay    rate  was   assumed   to   be  constant    in   the   entire 

stretch  of  the   river   under   study.      The   temperature  correction 

coefficient  was    1.106   for   the   scale-up  of   the  decay   rate   to   the 

design   temperature   conditions. 

4.3.3         Modell ing   Results: 

Table    3   shows   a    surmiary  of  the  modelling    results    for   the   24  options 

stated   above.      These   results   can   also   he   used   to  examine   the 

sensitivity   of  the  model    to   changes    in   various   design   parameters. 

The   input   parameters   Q   ,    Q,    temperature   and   pH  are   listed   in 

columns   2  to   5.      The    instream  average  concentration   at   the 

downstream  boundary   (i.e.   just   upstream  of  the   Kitchener   WPCP 

outfall)    is   given    in   column   6.      The   downstream  distance   along   the 

shoreline   x        ,    at  which  the  concentration   of  total    amnonia 
see 

becomes   c  ^   when   c     =   12.5  mg/L   (as    N),    is   listed    in   column   7. 
St  e  a/      \  /  ' 

The   next  column  shows   the   lateral    boundary  of  the   LUZ     q. ,    as   a 

fraction   of   the    river  discharge,    and   the   last   two   columns    include 

the   values   of  allowable   effluent   concentration,   c   .,    and   the 

shoreline  distance,   x        .      at   which   the   criterion    is  met. 

sceA 

(Note:   In  Table  3,  x    =  -999.0  indicates  that  x    exceeds 
^  see  see 

the  limiting  distance  of  20,487  m). 

The  following  general  features,  which  indicate  sensitivity  aspects 
as  well  as  the  relative  effects  of  various  management  options,  are 
evident  from  the  results  presented  in  Table  3: 


1.  When  c  =  12.5  mg/L  (as  N)  ,  x    decreases  due  to  a  decrease 
in  temperature,  pH  and  Q  ,  and  an  increase  in  Q. 
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TABLE  3  :   AMMONIA  IN  THE  GRAND  RIVER  -  SUMMARY  OF  PREDICTIONS 

rUHMARY  OF  RUNS  FOR  HANACEMENT  OPTIONS' 

AHHONrA  IN  MIXING  ZONE-r.RnND  RIVER  IiELQU  UATEPLOO  UPCP  • 

RUNÎ  OEFL  OPUF  TEMP     PH     CAUP      XSC£            OY/OT     CEA      XTCEA 

1.  0.307  4.24  15.0 

2.  0.3G7  4.24  15.0 

3.  0.307  4.24  20.0 

4.  0.307  4.24  20.0 

5.  G.3©7  4.24  25.0 

6.  0.307  4.24  2r.O 

7.  e.307  5.49  15.0 
■e.  0.307  8.49  15.0 
9.  9.307  8.49  20.0 

10.  0.307  8.49  20.0 

11.  e.307  8.49  25.0 

12.  e.307  8.49  25.0 

13.  O.Aee  4.24  15.0 

14.  0.6e0  4.24  15.0 

15.  ^.680  4.24  20.0 

16.  0.600  4.24  20.0 

17.  0.6ee  4.24  25.0 

18.  0.600  4.24  25.0 

19.  0.600  8.49  15. G 

20.  Ô.60G  8.49  15.0 

21 .  0.600  8.49  20.0 

22.  0.600  8.49  20.0 

23.  0.600  6.49  25. '.1 

24.  O.oGO  9.49  25.'.' 


7. 

8 

0. 

358 

5410. 

7 

0.20 
0.30 

0.40 

7.56 
12.17 
16.91 

2970.8 
5037.8 
7892.2 

8. 

3 

0. 

358 

1  941  7. 

0 

0.20 
0.30 
0.40 

2.48 
3.99 
3.54 

3054.7 
5036.3 
7092.2 

7. 

3 

0. 

205 

6855. 

1 

0.20 
0.30 
I.'.  40 

5  .  40 

9.  04 

12.83 

2789.6 
4725. 1 
7040.2 

8. 

3 

n . 

205 

17219. 

4 

0.20 
0.3G 
0.40 

1  .80 
3.01 
•1.27 

2738.3 
4726.3 
7040.2 

7  . 

9 

0. 

091 

6989. 

8 

0.20 
0.30 
0.40 

i.8f 

7.05 

1Û.3Ô 

2  :24.2 
4234. 1 
5975.0 

a. 

,3 

0. 

081 

14391 . 

8 

0.20 
0.30 
0.40 

1  .31 
2.40 
3.50 

2462. 1 
4231 .3 
5975.1 

7. 

.8 

0. 

231 

2579  , 

,0 

0.20 
0.30 
0.40 

14.45 
22.93 
31  .61 

3143.1 
5576.7 
8539.1 

3. 

3 

0. 

231 

1 1203. 

,3 

0.20 
0.30 
0.40 

4.74 

7.52 

10.36 

3255.1 
5575.9 
8539.1 

7. 

.8 

0. 

152 

3398. 

,0 

0.20 
0.30 
0.40 

10.27 
16.81 
23.52 

2975.5 
5040.4 

7722.0 

3. 

,3 

0. 

152 

12409, 

,3 

0.20 
0.30 
0.40 

3.42 
5.60 
7.83 

3064.9 
5038.5 
7722.1 

7, 

.8 

0. 

076 

4496. 

,  1 

0.20 
0.30 

0.40 

7.50 
12.97 
18.50 

2784.2 
4710,3 
6813.6 

8. 

,3 

0. 

076 

1 1843, 

,9 

0.20 
0.30 
0.40 

2.55 
4.41 
6.2" 

2733.3 
4711  .7 
6813.6 

7, 

.8 

0. 

,674 

1 1 621 , 

.9 

0.20 
0.30 

o.4e 

4.11 
6.61 
9.19 

3088.2 
5062.0 
7=64.7 

8. 

.3 

0. 

.674 

-999, 

.0 

0.20 
0.30 
0.40 

1  .35 

2.17 
3.01 

2967.0 
5064.2 
7964.7 

7 

.8 

0, 

.391 

12154 

n 

0.20 
0.30 
0.40 

2.94 
4.91 
6.95 

2774.5 
4756.6 
7106.5 

8 

.3 

0 

.391 

-999 

.0 

0.20 
0.30 
0.40 

0.98 
1.63 
2.31 

2804.3 
4755.9 
7106.5 

7 

.8 

0 

.  159 

1  1042 

.3 

0.20 
0.30 
0.40 

2.09 
3.78 
5.57 

2395.0 
4222.3 
6034.3 

8 

.3 

0 

.  15° 

19260 

.5 

0.20 
0.30 
0.40 

0.71 
1  .29 
1  .69 

2470.5 
4216.9 
6034.5 

-» 

.8 

0 

.  441 

5758 

.  6 

0.20 
0.30 
0.40 

7.64 
12.11 
16.69 

3280.8 
55=2.1 
9570.9 

8 

.  3 

0 

.  441 

-«9V 

.0 

0.20 
0.30 
0.40 

2.50 
3.9^ 
5.47 

32^1 .6 
5592.2 
8570.9 

7 

.8 

0 

.293 

7818 

.8 

0.20 
0.30 
0.40 

5.42 

8.07 

12.42 

2961 .9 
5050.1 
7763.0 

8 

.3 

0 

.2"3 

-9WO 

.0 

0.20 
0.30 
0.40 

1  .91 

2  .  95 
4.13 

2962.3 
5056.0 
7763.0 

.8 

0 

.  1  48 

9539 

.7 

0.20 
0.30 
0.40 

3.96 
6.85 
9.  75 

27oi .0 
4735. 7 
68  49..: 

8 

.3 

0 

.  1  46 

1  7=1  5 

.6 

0.20 
0.30 
0.4G 

1  .35 
2.33 
3.31 

2719.5 
4734.9 
6949.6 

2.  The  concentration  of  total  amionia  at  the  downstream  boundary 
(i.e.,  just  upstream  of  the  Kitchener  outfall)  decreases  due  to 
increases  in  temperature  and  Q,  and  a  decrease  in  Q  . 

3.  An  increase  in  the  lateral  boundary  of  the  LUZ  alone  results  in 
increases  in  both  c  .  and  x   .. 

4.  A  sole  increase  in  pH  would  decrease  c  .;  the  values  of 

eA 

X   .  are  not  affected  by  pH.  (Note:  Some  values  of  x   . 
sceA  ^   ^    ^  sceA 

are  seen  to  differ  due  to  a  change  in  pH  only;  this  is  due  to  a 
5%  difference  permitted  in  an  iterative  procedure). 

5.  As  temperature  increases,  both  c  .  and  x   .  decrease. 

6.  An  increase  in  the  nver  discharge,  Q,  would  result  in  increased 
Sa  ^"^  Seek- 


As  Q  increases,  c  .  decreases  and  x   .  increases. 
e  eA  sceA 


(Note:   In  some  cases,  the  values  of  x   »  are  seen  to 
decrease  with  an  increase  in  Q  ;  this  may  be  due  to  a  ! 
difference  permitted  in  the  iterative  procedure.) 


It  should  be  noted  that  the  foregoing  salient  features  are  of  a 
general  nature  valid  for  the  case  of  a  nonconservati ve  material 
following  first  order  decay  kinetics  (viz.,  chlorine,  phenol,  etc.); 
however,  the  features  involving  pH  effects  are  specific  to  the  case 
of  ammonia  only. 
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4.4 


Residual    Chlorine    in   the   Boyne   River   below  A1 1 i ston 


4.4.1         General    Description  of   Studies: 

Chlorine   is   widely  used   as   a  disinfecting   agent    in   wastewater 
treatment   practice    in   order   to  destroy  microorganisms   of   public 
health   concern.      In   recent  years,   the   discharge  of  chlorinated 
effluents    into    streams,    rivers   and    lakes   has   caused   concern   because 
of  the   potential    toxicity   of  chlorine   to   aquatic   biota.      The   MOE 
Water   Quality   Objective   for   total    residual    chlorine   (TRC)    is   0.002 
mg/L   (10). 

The   location   selected    for   this   study   is   a   2.9  km  segment   of  the 
Boyne  River   located   between   the  Alii ston   WPCP  outfall    and   the 
confluence   of  the   Boyne  with   the   Nottawasaga   River.      Effluent    from  a 
food   processing    industry   (Salada   Foods)    is  discharged    into   the   creek 
at   about    945  m  below  the   Al  I  i  ston   WPCP  outfall.      A  schematic    layout 
of   the    study    stretch  of   the    river    is    shown    in   Fig.    10. 

During   the   field    investigations,   data   on   channel    geometry,    flow 
rates,    time  of  travel    and   effluent  dispersion  characteristics  were 
collected.      Instream  residual    chlorine   levels  at   various   distances 
below  the  outfall    were  measured  during  a   preliminary   run.      Based   on 
these  data,   several    transects  were   established    in   the   stream  for 
sampling   purposes.      The   details  of   field  data   at   each   location, 
including   the   number   of   sampling    runs   and   parameters   measured,    are 
presented  by  Wisz,    Ellis   and    Inniss    (15).      Some   results  of  data 
analysis   and   modelling   studies   are   provided    in    previous   publications 
(4,    6,    8,    12).      A  brief   sumnary   of   the    field  data    is   presented    below. 


An    intensive   field    survey  was   carried   out    in   the   Boyne   River   during 
August    11-12,    1976.      During   the   survey,    the    st^eamflow    in    the   Boyne 
River   upstream  of  the   Alliston   WPCP  was  measured   at   0.78  m   /sec 
and   the   average   instream  temperature  was    21.5   C;   the   average 
effluent   flow    rate  was    39.6  L/sec.      The  average   concentrations   of 
total    residual    chlorine   and   chloride    ion   at   various   sampling   points 
for   the  duration   of   the   survey   are   suimarized    in   Table   4. 
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METERS 


FIGURE  10.  SCHEMATIC  LAYOUT  OF  STUDY  AREA  -  BOYNE  RIVER 
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TABLE   4:      SUr^ARY  OF   BOYNE   RIVER   FIELD   SURVEY  DATA 
AUGUST    11-12,    1976. 


Transect 

Distance 

Transverse 

Aver-age 

and 

below 

Distance 

below 

Avg^ 

Chlorides'^ 

Sampl ing 

Outfall 

Outfall 

Bank 

TRC 

as  CL- 

Station 

m 

m 

ug/L 

mg/L 

Control 

0 

13.7 

Effluent 

0 

1170 

131.0 

Transect   A 

21 

8.84b 

lOOC 

A  1/2 

0.381 

220 

47.0 

Al 

0.76 

18 

27.0 

A2 

1.52 

0 

16.2 

A3 

2.29 

0 

14.4 

A4 

4.88 

0 

13.5 

A5 

6.71 

0 

13.6 

Transect   B 

61 

10.70b 

16C 

Bl 

1.52 

23 

24.3 

B2 

3.05 

10 

20.4 

B3 

4.57 

2 

17.0 

84 

7.62 

0 

13.6 

B5 

9.45 

0 

13.4 

Transect   C 

183 

6.40b 

3C 

CI 

1.22 

6 

20.0 

C2 

2.44 

3 

19.0 

C3 

3.66 

1.4 

17.0 

C4 

4.88 

0.6 

15.0 

Salada   Foods 

945 

0 

21.0 

Outfall 

Transect   D 

1148 

0 

18.0 

Transect   E 

2897 

0 

17.0 

^       Average  concentration   of   five  samples   taken   over   a  duration   of  27 
hours    (except   three   samples   at   1   1/2). 

b       Channel    width   at   the   transect. 

*-       Mean   concentration    in    the   effluent   plume. 
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4.4.2         Data   Analysis   and    Parameter   Estimation: 

The  details  of   estimation   of  various   parameters   are  described    in 
previous   publications    (4,    6,    8).      The   salient   results  of  data 
analysis   are   sumnarized   below. 

The   simple  method   described   in    Section    3.3  has   been   used   to   estimate 
the   coefficient,   6    .      Using   the  chlonde    ion  data   presented    in 
Table   4,   the   net  maximum  concentration   at   each   transect  was 
determined    from  the   observed  maximum  concentration,  minus   the 

background   value;   the   value   of  c     is   6.33  mg/L  from  the  data   given 

2     ^ 
above.     The   plot   of   (c  /c_     )      versus    (x/b)    is   shown   in 

d     max 

Fig.    11,    from  which   the  value  of    B  is   found   to   be   0.0038,    based   on 
an   estimated   best-fit   line. 

The   available  data   on   cross-sectional    profiles   of  depth,   velocity 
and  concentration  of   TRC  were   insufficient   to   compute   the  mass   flux 
of  TRC   at   each  transect    (required   to  determine  decay   rates   from  a 
semi -1 ogarithmic  plot).      Hence,    an   alternative  method    involving 
determination   of  gross   decay   rates    (due   to   dispersion,   chemical 
reaction,    etc.)    and   apparent   attenuation    rates   (due   to  dispersion 
only),   has   been   used   to   calculate   the  decay   rates   of  TRC;   this 
method    is   described   below. 

The   gross   TRC   decay   rate  coefficient,   k    ,   due  to  mixing   and   all 
other  mechanisms,  was   determined    from  a    semi-log   plot   of   observed 
average   TRC   concentration    in    the   plume  vs.   distance,   shown   in 
Fig.    12.      An   examination   of   the  plot    indicates    that   the    rate  varies 
significantly   in    two  successive   reaches   as    indicated   by  the   slopes 
of  the   solid    lines.      The   estimated    rate   coefficients,    k    ,    for   the 
two   reaches   AB   and   BC   are   0.0104   and   0.0031   per   second    (to   base   e) , 
respectively,    at   the   observed  mean   temperature   of   21.5   C. 

The  width   of   the   effluent   plume   at   each   transect  was   estimated    from 
the   chloride   ion   distributions.      The   partial    discharge  within   the 
plume  width    at   each   transect  was    then    calculated   on   the   assumption 
that   the   discharge   per   unit   width   of  channel    is   the   same.      Treating 
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the   TRC   as   a  conservative  material,   the  average  TRC  concentration 

values  within    the   effluent   plume,   due   to  mixing   effects   only   and 

computed    from  mass   balance,   were;      177.69,    112.65,    70.38  and 

56.98  ug/L   at   Transects   A,    B,    C   and   D,    respectively.      On   the 

assumption   that   the   decrease    in   concentration   due   to  mixing    follows 

an   exponential    law,    the   apparen';   attenuation   rate   coefficients  due 

to  mixing    (k    )   are   found   to   be   0.0026  and   0.00088  per   second    (to 

base  e)    for   the   reaches   AB   and   BC,    respectively   (see  Fig.    12), 

representing   approximately  25%  and   28%  of  the  gross  decay   rate 

values.     The   corresponding  decay   rates   of  TRC    (k.)   attributable   to 

chemical    reactions   and   other  mechanisms,   determined    from  k,   =   (k.-k    ), 

d        ^   t     p^  ' 

are   0.0078  and   0.0022   per   second    (to   base   e)    at   21.5  C. 

4.4.3         Model    Validation   Study: 

Due  to   the   limitations   of  data   available,   the  model    validation   was 
carried  out    in   one    step   only    (instead   of   the   calibration   and 
verification   procedures).      Using   the   parameters   obtained    from  the 
survey  data   as    input   to   the  MIXCALBN  program     concentration 
predictions   were   obtained  with      6  ■    =  0.0038;   however,   the 
predictions  were   not    in   agreement  with   the   observations.      Then, 
several    runs   were   carried   out    by   selecting   different   combinations 
of    S    and   k     for   each  transect;   and   satisfactory  agreements 
between   the   observations   and   the  predictions,    shown    in   Fig.    13,   were 
obtained  with   the   following  combinations: 

Transect  s  ^d 


A  0.0005  0.007/sec. 

B  0.0025  0.007/sec. 

C  0.0025  0.0022/sec. 


The   predictions   shown    in   Fig.    13   are   based    on   the   approximation    that 
the   discharge   per   unit   width   of   channel    is    the   same  at   each 
transect.      Even  with   this    approximation,    the   predictions   of   the 
stream  tube   model    are    in   good   agreement   with   the   observations.      In 
contrast,    results   presented    in   previous   publications    (4,    6)    show 
that   the   predictions   of  the  models   based   on   width   as   the   lateral 
variaole  differ   considerably   from   the   observations   at    some   transects. 
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4.4.4       Critical    Point   Method   Predictions: 

Using   the   parameters   determined    in    the  model    validation    step,    the 
transverse   distributions   of   TRC  were   predicted   at    several    closely 
spaced  distances   below  the  outfall.      From  these  predictions,    the 
longitudinal    distributions   of   TRC   along   q.     =  0.;.'5Q  and   q,     =  0.35Q 
were  plotted   as    shown    in   Fig.    14.      The   occurrence   of   the   critical 
point    is   clearly  demonstrated    in   these   plots.      The  critical    concen- 
trations of   TRC   are   seen   to   be   higher   than   the   PWQO  of   2  ug/L.      For 
comparison   purposes,   the   values   of   TRC   at   the   appropriate   values   of 
q.    interpolated    from  the   observed   cross-   sectional    distributions 
at   Transects   B   and  C   are   also   shown   on   the  plots. 

For  the   same  design   conditions   as   above   -  viz.,    c     =   1170  ug/L, 
c^    =  2  ug/L,    Q^   =  0.04  m'^/s,   and   Q  =  0.82  m'^/s,   the  critical 
co-ordinates   and   allowable  effluent   TRC   concentrations  were 
calculated   by  the   graphical    method  described   in   Section   2.5,    using 

estimated   average  values    r     =  0.0025,   k       =  0.007  per   second   at 

^  av 

21.5°C   (to   base   e),   B   =  9.9  m,    and   U   =  0.226  m/s.      The  completely 
mixed   instream  average  concentration,   c     =  57.07   ug/L.     The 
results   for   p.    (i.e.,    q./Q)    =  0.25  and   p.    =  0.35  are   given   in 
Table   5.      For   comparison   purposes,   the   values   of  x,  ,   c.    and 
c    .   predicted    from  the   MIXAPPLN   program  are   also   given    in    Table   5 
(inside   parentheses).      The  predictions  of  the  two  methods   are   seen 
to   be    in   reasonable   agreement. 

TABLE  5.      RESULTS   OF   CRITICAL    POINT  METHOD 


Pl 

^L 

\ 

V 

\ 

^L 

'^eA 

X   reduction 

= 

m 

ug/L 

ug/L 

100(c^-c^,). 

's 

0.25 

122. 

7 

1.2 

0.30 

38.7 
(35.0) 

17.1 
(19.5) 

136.8 
(120.0) 

88.3 
(89.7) 

0.35 

122. 

7 

1.7 

0.088 

54.9 
(57.0) 

5.0 
(6.0) 

468.0 
(390.0) 

60.0 

(56.7) 
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With  the  design  parameters  listed  in  the  previous  paragraph, 
computations  based  on  the  instantaneous  compl ete  mixing  (ICM) 
assumption  would  result  in  an  allowable  effluent  TRC  concentration 
of  41  ug/L,  requiring  a  reduction  of  96.5%.   In  contrast,  the  design 
based  on  the  LUZ  concept  with  p.  =  0.35  requires  a  reduction  of 
about  60%  only.  The  results  of  examples  presented  in  another 
publication  (7)  indicate  that  in  some  cases,  the  ICM  assumption  or 
the  dilution  ratio  concept  may  lead  to  a  LUZ  with  too  large  a 
lateral  boundary  (and  consequently,  too  small  a  zone  of  passage)  due 
to  an  underestimation  of  treatment  required  in  comparison  to  the 
predictions  based  on  the  LUZ  concept;  whereas,  in  some  other 
situations,  the  predictions  based  on  the  former  method  may  result  in 
treatment  requirements  that  are  too  stringent  in  comparison  to  those 
of  the  latter  method.  Thus,  in  the  case  of  discharge  from  a  bank 
outfall,  it  is  inappropriate  to  determine  the  allowable  effluent 
concentration  of  a  pollutant  based  on  the  instantaneous  complete 
mixing  assumption. 
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V.      SUftlARY  AND   CONCLUSIONS 

This   report   deals  with  an   assessment   procedure   for   prediction   of 
water  quality   in  mixing  zones  of   shallow   rivers   below  bank 
outfalls.      It    is  mainly  an   abridged   version  of  the   MOE   publication, 
Water   Resources   Paper  #14   (Ref.   8). 

Descriptions   of  modelling   concepts    including   expressions    for   the 
critical    point,    and   allowable   effluent   concentration   of  a   pollutant 
to  meet   a  given    instream  objective  are   presented.      The  computations 
can   be   performed  with   the   aid  of  charts   presented    in   this    report   and 
desk-top  calculations,   or   through  the   use  of  computer  programs 
documented    in  Water  Resources   Paper  #14   (Ref.   8).      The   field   survey, 
data  analysis   and   parameter  estimation   procedures   are  described.      A 
step-by-step  design   procedure    is  outlined   to   perform  computations 
with   the   aid  of  charts   and  desk-top  calculators,   and    is   illustrated 
by   an   example.      The   application  of  a   computer  model    to  evaluate   the 
impact   of  typical    management   options   on   the  mixing   zone  boundaries 
and   allowaole   effluent   concentrations    is   also   illustrated   through 
examples. 

The    following   are   the   general    conclusions   of  this   study: 

1.  It    is  more   appropriate   to  designate   the   lateral    boundary  of   LUZ 
and   ZOP   as    fractions   of  total    discharge   in   river   instead   of 
channel   widths,    particularly    in   shallow   rivers  of  nonuniform 
cross-sectional    geometry. 

2.  The   steady  state   stream  tube  model,   based   on  modified   analytical 
solutions   of   the   2-D   convection-dispersion   equation,    can   be 
applied   to   natural    streams   and    rivers    in   which   tracer   and 
pollutant   concentration   profiles   are   continuous   at   successive 
transects. 

3.  An    increase    in   the    lateral    boundary   of  a    limited   use   zone 
results    in    an    increase    in   the   allowable   effluent   concentration; 
however,   the   associated    longitudinal    spread   along   the   discharge 
shoreline  would   also    increase. 


-   56   - 


4.  When  the  discharge  in  a  river  is  increased,  the  allowable 
effluent  concentration  as  well  as  the  maximum  longitudinal 
spread  would  increase. 

5.  As  temperature  increases,  the  allowable  effluent  concentration 
and  the  maximum  longitudinal  spread  decrease. 

6.  In  the  case  of  discharge  from  a  bank  outfall,  it  is 
inappropriate  to  determine  the  allowable  effluent  concentration 
of  a  pollutant  based  on  the  instantaneous  complete  mixing 
assumption  (as  well  as  the  dilution  ratio,  streamflow 
rate/effluent  flow  rate)  due  to  the  following  reasons:  (a)  in 
some  cases  it  may  lead  to  a  LUZ  with  too  large  a  lateral 
boundary  (and  consequently,  too  small  a  zone  of  passage)  due  to 
an  underestimation  of  treatment  required;  and  (b)  in  some  other 
cases,  it  may  result  in  too  small  a  LUZ  due  to  the  prediction  of 
too  stringent  treatment  requirements  in  comparison  to  those 
based  on  the  LUZ  concept. 
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